TTK-12-05, TUM-HEP-852/12 



Dark Matter, Baryogenesis and Neutrino Oscillations 
from Right Handed Neutrinos 

Laurent Canetti a , Marco Drewes b ' c , Tibor Frossard d , Mikhail Shaposhnikov a 

a ITP, EPFL, CH-1015 Lausanne, Switzerland 

b Institut fur Theoretische Teilchenphysik und Kosmologie, 

RWTH Aachen, D-52056 Aachen, Germany 

c Physik Department T31, Technische Universitat Miinchen, 

James Franck Strafie 1, D- 85748 G arching, Germany 

Q-I d Max- Planck- Institut fur Kernphysik, Saupfercheckweg 1, 69117 Heidelberg, Germany 

D 
43 



Abstract 

We show that, leaving aside accelerated cosmic expansion, all experimental data in high 
energy physics that are commonly agreed to require physics beyond the Standard Model can 
be explained when completing it by three right handed neutrinos that can be searched for 
using current day experimental techniques. The model that realizes this scenario is known as 
Neutrino Minimal Standard Model (VMSM). In this article we give a comprehensive summary 
of all known constraints in the z/MSM, along with a pedagogical introduction to the model. 
We present the first complete quantitative study of the parameter space of the model where 
no physics beyond the ^MSM is needed to simultaneously explain neutrino oscillations, dark 
matter and the baryon asymmetry of the universe. This requires to track the time evolution 
of left and right handed neutrino abundances from hot big bang initial conditions down to 
temperatures below the QCD scale. We find that the interplay of resonant amplifications, CP- 
violating flavor oscillations, scatterings and decays leads to a number of previously unknown 
constraints on the sterile neutrino properties. We furthermore re-analyze bounds from past 
collider experiments and big bang nucleosynthesis in the face of recent evidence for a non- 
zero neutrino mixing angle 613. We combine all our results with existing constraints on dark 
matter properties from astrophysics and cosmology. Our results provide a guideline for future 
experimental searches for sterile neutrinos. A summary of the constraints on sterile neutrino 
masses and mixings has appeared in [l[ . In this article we provide all details of our calculations 
and give constraints on other model parameters. 



Contents 

1 Introduction 

2 The i/MSM fjj 

2.1 Mass- and Flavor Eigcnstates 5 

2.2 Benchmark Scenarios 6 

2.3 Effective Theory of Lepton Number Generation 7 

2.4 Thermal History of the Universe in the ^MSM 

2.5 Parameterization 

2.6 "Fine Tunings" and the Constrained vMSM 

2.6.1 Baryogenesis 

2.6.2 Dark Matter Production 

3 Experimental Searches and Astrophysical Bounds 

3.1 Existing Bounds 

3.1.1 Seesaw Partners N2 and A3 

3.1.2 Dark Matter Candidate Ni 

3.2 Future Searches 

3.2.1 Dark Matter Candidate Ni 

3.2.2 Seesaw Partners iV2,3 

4 Kinetic Equations 

4.1 Short derivation of the Kinetic Equations 

4.2 Computation of the Rates 

4.2.1 Baryogenesis 

4.2.2 Dark Matter Production 

5 Baryogenesis from Sterile Neutrino Oscillations 

6 Late Time Lepton Asymmetry and Dark Matter Production 

7 DM, BAU and Neutrino Oscillations in the i/MSM 

8 Conclusions and Discussion 
A Kinetic Equations 

A.l How to characterize the Asymmetries 

A. 2 Effective Kinetic Equations 

A. 3 The Effective Hamiltonian 

A. 3.1 Dispersive Part H 

A. 3. 2 Dissipative Part 

A. 3. 3 The remaining Rates 

A. 4 Uncertainties 

B Connection to Pseudo-Dirac Base 

C How to characterize the lepton asymmetries 

D Low temperature Decay Rates for sterile Neutrinos 

D.l Semileptonic decay 

D.2 Leptonic decay 



2 



1 Introduction 



The Standard Model of particle physics (SM), together with the theory of general relativity 
(GR), allows to explain almost all phenomena observed in nature in terms of a small number 
of underlying principles - Poincare invariance, gauge invariance and quantum mechanics - and a 
handful of numbers. In the SM these are 19 free parameters that can be chosen as three masses for 
the charged leptons, six masses, three mixing angles and one CP violating phase for the quarks, 
three gauge couplings, two parameters in the scalar potential and the QCD vacuum angle. Three 
leptons, the neutrinos, remain massless in the SM and appear only with left handed chirality. 
GR adds another two parameters to the barcode of nature, the Planck mass and the cosmological 
constant. 

Despite its enormous success, we know for sure that the above is not a complete theory of 
nature for two reasons^- On one hand, it treats gravity as a classical background for the SM, 
which is a quantum field theory. Such description necessarily breaks down at energies near the 
Planck scale Mp and has to be replaced by a theory of quantum gravity. We do not address this 
problem here, which is of little relevance for current and near-future experiments. On the other 
hand, the SM fails to explain a number of experimental facts. These are neutrino oscillations, 
the observed baryon asymmetry of the universe (BAU), the observed dark matter (DM) and 
the accelerated expansion of the universe today. In addition there is a number of cosmological 
problems (e.g. flatness and horizon problem). These can be explained by cosmic inflation, another 
phase of accelerated expansion in the universe's very early history, for which the SM also cannot 
provide a mechanism. To date, these are the only confirmed empirical proofs of physics beyond 
the SrvH In this article we argue that, leaving aside accelerated cosmic expansion, all of them 
may be explained by adding three right handed (sterile) neutrinos to the SM that can be found 
in experiments. 

The model in which this possibility can be realized is known as Neutrino minimal Standard 
Model (i/MSM) @,3. The fMSM is an extension of the SM that aims to explain all experimental 
data with only minimal modifications. This in particular means that there is no modification 
of the gauge group, the number of fermion families remains unchanged and no new energy scale 
above the Fermi scale is introduced^). The matter content is, in comparison to the SM, comple- 
mented by three right handed counterparts to the observed neutrinos. These are singlet under 
all gauge interactions. Over the past years, different aspects of the i/MSM have been explored 
using cosmological, astrophysical and experimental constraints Moreover, it was 



suggested that cosmic inflation [3lM33l| and the current accelerated expansion [34H37I] may also 
be accommodated in this framework by modifications in the gravitational sector, which we will 
not discuss her^|. However, though the abundances of dark and baryonic matter have been es- 



x We do not address theoretical issues of "aesthetic" nature such as fine tuning in the context of the hierarchy 
problem, the strong CP problem and the flavor structure. They may be interpreted as hints for new physics, but 
could also simply represent nature's choice of parameters. 

2 We leave aside all experimental and observational anomalies that have not lead to a claim of detection of new 
physics, i.e. may be explained within the SM or by systematic errors. This includes the long standing problem of 
the muon magnetic moment, the inconclusive results of different direct DM searches as well as various anomalies 
of limited statistical significance. 

3 Because of this the technical hierarchy problem may be absent in the i/MSM because no new states with 
energies between the electroweak and the Planck scale are required 0, [B| . 

4 Inflation can be realized without modification of the gravitational interaction by adding an extra scalar to the 
i/MSM [3|| (see also [TBI. |39|. |40||). This inflaton can be light enough to be detected in direct searches. 
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timated individually in the framework of the i/MSM, to date it has not been verified that there 
is a range of right handed neutrino parameters for which they can be explained simultaneously, 
in particular for experimentally accessible sterile neutrinos. In this article we present detailed 
results of the first complete quantitative study to identify the range of parameters that allows 
to simultaneously explain neutrino oscillations, the observed DM density Odm and the observed 
BAU 4l|], responsible for today's remnant baryonic density £Ib- We in the following refer to 
this situation, in which no physics beyond the i/MSM is required to explain these phenomena, as 
scenario I. In this scenario DM is made of one of the right handed neutrinos, while the other 
two are responsible for baryogenesis and the generation of active neutrino masses. We also study 
systematically how the constraints relax if one allows the sterile neutrinos that compose DM to 
be produced by some mechanism beyond the i/MSM (scenario II). Finally, we briefly comment 
on a scenario III, in which the fMSM is a theory of baryogenesis and neutrino oscillations only, 
with no relation to DM. A more precise definition of these scenarios is given in section 12.21 Only 
scenarios I and II are studied in this article, which is devoted to the i/MSM as the common origin 
of DM, neutrino masses and the BAU. While scenario II has previously been studied in (22|, the 
constraints coming from the requirement to thermally produce the observed &DM i n scenario I 
are calculated for the first time in this work. We combine our results with bounds coming from 
big bang nucleosynthesis (BBN) and direct searches for sterile neutrinos, which we re-derived in 
the face of recent data from neutrino experiments (in particular 613 7^ 0). 

Centerpiece of our analysis is the study of all lepton numbers throughout the evolution of 
the early universe. As will be explained below, in the i/MSM lepton asymmetries are crucial 
for both, baryogenesis and DM production. We determine the time evolution of left and right 
handed neutrino abundances for a wide range of sterile neutrino parameters from hot big bang 
initial conditions at temperatures T ^> Tew ~ 200 GeV down to temperatures below the QCD 
scale by means of effective kinetic equations. They incorporate various effects, including thermal 
production of sterile neutrinos from the primordial plasma, coherent oscillations, back reaction, 
washouts, resonant amplifications, decoherence, finite temperature corrections to the neutrino 
properties and the change in effective number of degrees of freedom in the SM background. 
Many of these were only roughly estimated or completely neglected in previous studies. The 
various different time scales appearing in the problem make an analytic treatment or the use 
of a single CP-violating parameter impossible in most of the parameter space. Most of our 
results are obtained numerically. However, the parametric dependence on the experimentally 
relevant parameters (sterile neutrino masses and mixings) can be understood in a simple way. 
Furthermore, we discover a number of tuning conditions that can be understood analytically and 
allow to reduce the dimensionality of the parameter space. 

We find that there exists a considerable fraction of the i/MSM parameter space in which the 
model can simultaneously explain neutrino oscillations, dark matter and the baryon asymmetry of 
the universe. This includes a range of masses and couplings for which the right handed neutrinos 
can be found in laboratory experiments The main results of our study, constraints on sterile 
neutrino masses and mixings, have previously been presented in Q. In this article we give details 
of our calculation and constraints on other model parameters, which are not discussed in 

The remainder of this article is organized as follows. In Section [2] we overview the i/MSM, its 
parametrization, and describe the universe history in its framework, including baryogenesis and 
dark matter production. In Section [3] we discuss different experimental and cosmological bounds 
on the properties of right-handed neutrinos in the z^MSM. In Section 0] we formulate the kinetic 
equations which are used to follow the time evolution of sterile neutrinos and active neutrino 
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flavors in the early universe. In Section [5] we present our results on baryogenesis in scenario II. In 
Section [6] we study the generation of lepton asymmetries at late times, essential for thermal dark 
matter production in the z^MSM. In Section [7] we combine the constraints of the two previous 
Sections and define the region of parameters where scenario I can be realized, i.e. the z^MSM 
explains simultaneously neutrino masses and oscillations, dark matter, and baryon asymmetry of 
the universe. In Section[8]we present our conclusions. In a number of appendices we give technical 
details on kinetic equations (jAj), on the parametrization of the z^MSM Lagrangian (JB]), on different 
notations to describe lepton asymmetries ([C]) and on the decay rates of sterile neutrinos (iDj) . 

2 The i/MSM 

The z^MSM is described by the Lagrangian 

C-vMSM = £>SM + WR$V R - L L Fv R § - Vr-F^Ll^ 

- - (j^Mm vr + VrMI v c r ). ( 1 ) 

Here we have suppressed flavor and isospin indices. Csm is the Lagrangian of the SM. F is a 
matrix of Yukawa couplings and Mm a Majorana mass term for the right handed neutrinos ur. 
Ll = (vL,eL) T are the left handed lepton doublets in the SM and $ is the Higgs doublet. We 
chose a basis where the charged lepton Yukawa couplings and Mm are diagonal. The Lagrangian 
CD) is well-known in the context of the seesaw mechanism for neutrino masses [42| and leptogenesis 
[431 ] . While the eigenvalues of Mm in most models are related to an energy scale far above the 
electroweak scale, it is a defining assumption of the i/MSM that the observational data can be 
explained without involvement of any new scale above the Fermi one. 

2.1 Mass- and Flavor Eigenstates 

For temperatures T < M\y below the mass of the W-boson we can in good approximation replace 
the Higgs field $ by its vacuum expectation value v = 174 GeV. Then (pQ) can be written as 

C = C S m + wrjQvrj - (mn)aiW, a ^R,i ~ (m* D ) a iV^jv L>a 

-\{{M m )ij^^r,j + (M M )h^R,j) (2) 

with the Dirac mass matrix mr> = Fv. When the eigenvalues of Mm are much larger than those 
of m,Di the seesaw mechanism naturally leads to light active and heavy sterile neutrinos. This 
hierarchy is realized in the z^MSM. 

In vacuum there are two sets of mass eigenstates; on one hand active neutrinos Vj with masses 
mi, which are mainly mixings of the SU(2) charged fields vl, 

Pl^ = (ui((i- lee^j v L - ev%) y, (3) 

with a j = (mr)M M l ) a i, and on the other hand sterile neutrino^ Nj with masses Mj, which are 
mainly mixings of the singlet fields Ur, 

PrNj = (ul ( (i - \e T e*^ u R + e T uij ) . (4) 

5 In 0] the notation is slightly different and the letter "iVj" does not denote mass eigenstates. 
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Here Pr,l are chiral projectors and Nj (v,) are Majorana spinors, the left chiral (right chiral) part 
of which is fixed by the Majorana relations Nf = Nj and Vj = v?. The matrix Un diagonalises 
the sterile neutrino mass matrix Mn defined below. The entries of the matrix 8 determine the 
active-sterile mixing angles. 

The neutrino mass matrix can be block diagonalized. At leading order in the Yukawa couplings 
F one obtains the mass matrices 

m„ = -9M M 9 T , (5) 
M N = M M + ^9M M + M^e T e*) . (6) 

The mass matrices m v and Mn are not diagonal and lead to neutrino oscillations. While there is 
very little mixing between active and sterile flavors at all temperatures of interest, the oscillations 
between sterile neutrinos can be essential for the generation of a lepton asymmetry. m„ can be 
parameterized in the usual way by active neutrino masses, mixing angles and phases, m u = 
i/„diag(mi, m^, m^)U^ . In the basis where the charged lepton Yukawas are diagonal, U v is 
identical to the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) lepton mixing matrix. 

The physical sterile neutrino masses Mj are given by the eigenvalues of M n Mn- In the seesaw 
limit Mtv is almost diagonal and they are very close to the entries of Mm- We nevertheless need 
to keep terms O(0 2 ) because the masses M2 and M3 are degenerate in the i/MSM, see section 
12.61 and the mixing of the sterile neutrinos amongst each other may be large despite the 
seesaw-hierarchjU. This mixing is given by the matrix Un, which can be seen as analogue to U v . 
It is worth noting that due to ([6]) the matrix Un is real at this order in F. The experimentally 
relevant coupling between active and sterile species is given by the matrix 9 witrfl 

Qai = (eU N ) aI = (mnM^UN)^. (7) 

In practice, experiments to date cannot distinguish the sterile flavors and are only sensitive to 
the quantities 

ul = J2 @ -ie* a i = ^20ai0* a i. (8) 

/ I 

Therefore Un, and hence the sterile-sterile mixing and the coupling of individual sterile flavors 
to the SM, cannot be probed in direct searches. 

2.2 Benchmark Scenarios 

The notation introduced above allows to define the scenarios I-III introduced in the introduction 
more precisely. 

• In scenario I no physics beyond the z^MSM is needed to explain the observed Qdai, neu- 
trino masses and £Ib- DM is composed of thermally produced sterile neutrinos N±. N2 
and A3 generate active neutrino masses via the seesaw mechanism, and their CP-violating 
oscillations produce lepton asymmetries in the early universe. The effect of N\ on neutrino 
masses and lepton asymmetry generation is negligible because its Yukawa couplings F a \ are 

6 It turns out that the region where Un is close to identity phenomenologically is the most interesting, see section 

m 

7 The fact that matrix appearing in Q is U^ N 9 T = (Q*Un)^ rather than = U N 9^ is due to the fact that the 
Ni couple to vl,c, but overlap with v c L a . 
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constrained to be tiny by the requirement to be a viable DM candidate, c.f. section 13.1.21 
The lepton asymmetries produced by iVjj 3 are crucial on two occasions in the history of 
the universe: On one hand the asymmetries generated at early times (T > 140 GeV) are 
responsible for the generation of a BAU via flavored leptogenesis, on the other hand the late 
time asymmetries (T ~ 100 MeV) strongly affect the rate of thermal N\ production. Due 
to the latter the requirement to produce the observed £Idm imposes indirect constraints on 
the particles iV^. There are determined in sections [6] and [7] and form the main result of 
our study. 

• In scenario II the roles of iV^ and N\ are the same as in scenario I, but we assume that 
DM was produced by some unknown mechanism beyond the z^MSM. The astrophysical 
constraints on the N\ mass and coupling equal those in scenario I. are again required 
to generate the active neutrino masses via the seesaw mechanism and to produce sufficient 
flavored lepton asymmetries at T ~ 140 MeV to explain the BAU. However, there is no need 
for a large late time asymmetry. This considerably relaxes the bounds on ./V23. Scenario II 
is studied in detail in section [5j 

• In scenario III the z^MSM is not required to explain DM, i.e. it is considered to be a theory 
of neutrino masses and low energy leptogenesis only. Then all three Nj can participate in 
the generation of lepton asymmetries. This makes the parameter space for baryogenesis 
considerably bigger than in scenarios I and II, including new sources of CP violation. We 



do not study scenario III in this work, some aspects are discussed in 30]. 



2.3 Effective Theory of Lepton Number Generation 

In scenarios I and II the lightest sterile neutrino N\ is a DM candidate. In this article we focus on 
those two scenarios. If N\ is required to compose all observed DM, its mass M\ and mixing are 
constrained by observational data, see section [3l Its mixing is so small that its effect on the active 
neutrino masses is negligible. Note that this implies that one active neutrino is much lighter than 
the others (with mass smaller than O(10" 5 ) eV Q). Finding three massive active neutrinos with 
degenerate spectrum would exclude the z^MSM with three sterile neutrinos as common and only 
origin of active neutrino oscillations, dark matter and baryogenesis. N\ also does not contribute 
significantly to the production of a lepton asymmetry at any time. This process can therefore 
be described in an effective theory with only two sterile flavors A^. In the following we will 
almost exclusively work in this framework. To simplify the notation, we will use the symbols 
Mtv and Un for both, the full (3 x 3) mass matrix and mixing matrices defined above and the 
(2x2 and 3x2) sub-matrices that only involve the sterile flavors 7 = 2,3, which appear in the 
effective theory. The mixing between N\ and ^2,3 is negligible due to the smallness of F a %, which 
is enforced by the seesaw relation ([5]) and the observational bounds on M\ summarized in Section 
13.1.21 The effective ^2,3 mass matrix can be written as 

M N = Ml 2X 2 + AM ( j3 + M~ 1 Re(m| ) m D ), (9) 

where 03 is the third Pauli matrix and we chose the parameterization Mm = diag(M — AM, M + 
AM). This equality holds because we chose Mm real and diagonal. The physical masses M2 and 
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zero abundance of N 



thermal production of N 2 , N 3 
— lepton asymmetry generated 



200 Gev Tew Electroweak Symmetry Breaking 

— ► lepton Asymmetry converted to baryon asymmetry 

T + N 2 , N 3 reach equilibrium 

— ► lepton asymmetry washed out 

fewGeV T. N 2 , N 3 freeze out 

— • lepton asymmetry generated 



Id N 2 , N 3 decay 

— lepton asymmetry generated 

100 Mev Tdm resonant Ni Dark Matter production 



Figure 1: The thermal history of the universe in the vMSM. 



M3 are given by the eigenvalues of M^. They read 



M 2i3 
M 

(5M) 2 



M±5M 
1 

2M" 

( Re K m °/33 



M + ^rrRe ( tr ( m) D m£> 



\ 2 1 

Re (m ] D m D ^ J + AM J + Re (m) D m D 
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(10) 

(11) 

(12) 



For all parameter choices we are interested in M ~ M holds in very good approximation. The 
masses M2 ; 3 are too big to be sensitive to loop corrections. In contrast, the splitting 5M can be 
considerably smaller than the size of radiative corrections to M<i$ 44] • The above expressions 
have a different shape than those given in 0] because we use a different base in flavor space, see 
appendix iBl 

These above formulae hold for the (zero temperature) masses in the microscopic theory. At 
finite temperature the system is described by a thermodynamical ensemble, the properties of 
which can usually be described in terms of quasiparticles with temperature dependent dispersion 
relations. We approximate these by temperature dependent "thermal masses". 

2.4 Thermal History of the Universe in the z^MSM 

Apart from the very weakly coupled sterile neutrinos, the matter content of the i/MSM is the same 
as that of the SM. Therefore the thermal history of the universe during the radiation dominated 
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era is similar in both models. Here we only point out the differences that arise due to the presence 
of the fields vr, see figure [H They couple to the SM only via the Yukawa matrices F, which are 
constrained by the seesaw relation. For sterile neutrino masses below the electroweak scale, the 
abundances are too small to affect the entropy during the radiation dominated era significantly. 
However, the additional sources of CP-violation contained in them have a huge effect on the 
lepton chemical potentials in the plasma. 



Baryogenesis The z^MSM adds no new degrees of freedom to the SM above the electroweak 
scale. As a consequence of the smallness of the Yukawa couplings F, the Nj are produced 
only in negligible amounts during reheating Therefore the thermal history for T S> Tew 

closely resembles that in the SMu- The sterile neutrinos have to be produced thermally from the 
primordial plasma in the radiation dominated epoch. During this non-equilibirum process, all 



Sakharov conditions 45(] are fulfilled: Baryon number is violated by SM sphalerons [461 ] . and 



the oscillations amongst the sterile neutrinos violate CP [47]]. Source of this CP-violation are 
the complex phases in the Yukawa couplings F a j. Due to the Majorana mass Mm neither the 
individual (active) leptonic currents, defined in ([93]) and ([94"|) . nor the total lepton number are 
strictly conserved. However, for T>M the effect of the Majorana masses is negligible. Though 
the neutrinos are Majorana particles, one can define neutrinos and antineutrinos as the two 
helicity states, transitions between which are suppressed at T > M. We will in the following 
always use the terms "neutrinos" and "antineutrinos" in this sense. 

In scenarios I and II the abundance of iVi remains negligible until T ~ 100 MeV because of 
the smallness of its coupling that is required to be in accord with astrophysical bounds on DM, 
see Section [3.1.21 on the other hand, are produced efficiently in the early universe. During 



this process flavored "lepton asymmetries" can be generated [47|. reach equilibrium at a 

temperature T + 0]. Though the total lepton number (|93p at T+ 3> M is very small, there are 
asymmetries in the above helicity sense in the individual active and sterile flavors. Sphalerons, 
which only couple to the left chiral fields, can convert them into a baryon asymmetry. The 
washout of lepton asymmetries becomes efficient at T < T + . It is a necessary condition for 
baryogenesis that this washout has not erased all asymmetries at Tew, which is fulfilled for 
TV > Tew- The BAU at T ~ Tew can be estimated by today's baryon to photon ratio, see 



41[ for a recent review. A precise value can be obtained by combining data from the cosmic 



microwave background and large scale structure [48[] 



1]B = (6.160 ±0.148) • 10" 10 . (13) 

The parameter t/b is related to the remnant density of baryons Qb, in units of the critical 
density, by £Ib — ??b/(2.739 • 10~ 8 h 2 ), where h parameterizes today's Hubble rate Hq = WOh 
(km/s)/Mpc. In order to generate this asymmetry, the effective (thermal) masses M2(T) and 
Mz(T) of the sterile neutrinos in the plasma need to be quasi-degenerate at T > Tew, see section 



After N2 and A3 reach equilibrium, the lepton asymmetries are washed out. This washout 
takes longer than the kinetic equilibration, but it was been estimated in 0] that no asymmetries 



survive until A2 i 3-freezeout at T = T_. In [49j] it has been suggested that some asymmetry may 



If a non-minimal coupling of the Higgs field $ to gravity is introduced in the z^MSM, $ can play the role of the 
inflaton. Though this way no fields are added, the thermal history at very early times (during reheating) changes 
due to a non-minimal coupling to curvature, see [33] ■ Here we assume an initial state without Ni at T > Tew- 
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be protected from this washout by the chiral anomaly, which transfers them into magnetic fields. 
Here we take the most conservative approach and assume that no asymmetry survives between 
T + and T_. Around T = T_, the interactions that keep iV^ in equilibrium become inefficient. 
During the resulting freezeout the Sakharov conditions are again fulfilled and a new asymmetries 
are generated. Even later, a final contribution to the lepton asymmetries are added when the 
unstable particles decay at a temperature T^. 



DM production The abundance of the third sterile neutrino N\ in scenario I remains below 
equilibrium at all times due to its small Yukawa coupling. In absence of chemical potentials, 
the thermal production of these particles (Dodelson-Widrow mechanism j5o| ) is not sufficient to 
explain all dark matter as relic Ai abundance if the observational bounds summarized in section 
[3] are taken into consideration. However, in the presence of a lepton asymmetry in the primordial 
plasma, the dispersion relations of active and sterile neutrinos are modified by the Mikheyev- 
Smirnov-Wolfenstein effect (MSW effect) [5l| . The thermal mass of the active neutrinos can 
be large enough to cause a level crossing between the dispersion relation for active and sterile 
flavors at Tdm, resulting in a resonantly enhanced production of N\ [3] (resonant or Shi-Fuller 



mechanism |52j]). This mechanism requires a lepton asymmetry \fJL a \ > 8 • 10 to be efficient 
enough to explain the entire observed dark matter density SI dm i n terms of N\ relic neutrinos 
[l9| | . Here we have characterized the asymmetry bjl^l 



Ma = — , (14) 
s 



where s is the entropy density of the universe and n a the total number density (particles minus 
antiparticles) of active (SM) leptons of flavor a. The relations between \i a defined in (j!4[) and 
other ways to characterize the asymmetry (e.g. the chemical potential) are given in appendix O 



Cosmological constraints Thus, in scenario I there are two cosmological requirements related 
to the lepton asymmetry that have to be fulfilled to produce the correct SIb and SI dm within the 
z/MSM: 

i) jjL a ~ 10~ 10 at Tew ~ 200 GeV for successful baryogenesis and 

ii) \/j, a \ > 8 • 10 -6 at Tdm f° r dark matter production. 

In scenarios I and II the asymmetry generation in both cases relies on a resonant amplification 
and quasi-degeneracy of M2 and M3 , which we discuss in section 12.61 This may be considered as 
fine tuning. On the other hand, the fact that the BAU (and thus the baryonic matter density SIb) 
and DM production in the SM both rely on essentially the same mechanism may be considered 
as a hint for an explanation for the apparent coincidence SIb ~ SIdm, though the connection is 
not obvious as SIb and Si dm also depend on other parameters. 

In scenario II only the condition i) applies. The resulting constraints on the A2,3 properties 
have been studied in detail in [22I ]. In section [5] we update this analysis in the face of recent data 
from neutrino experiments, in particular evidence for an active neutrino mixing angle 813 7^ 



53M55I] . In section [6] we include the second condition and study which additional constraints 



come from the requirement \fi a \ > 8 • 10 6 at Tdm- Previous estimates suggest Tdm ~ 100 



Note that fj, a is not a chemical potential, but an abundance (or yield). We chose the symbol fi for notational 
consistency with The relation of fi a to the lepton chemical potential (i.<x is given in appendix [C] 
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MeV < Tqcd and T_ < M w @, [Hj, where M w is the mass of the W^-boson and Tqcd the 
temperature at which quarks form hadrons. 

Though we are concerned with the conditions under which N\ can explain all observed dark 
matter, the N\ will not directly enter our analysis because the lepton asymmetry that is necessary 
for resonant N\ production in scenario I is created by ^2,3. Instead, we derive constraints on 
the properties of N 2j 3, which can be searched for in particle colliders. N\, in contrast, cannot be 
detected directly in the laboratory due to its small coupling. However, the N\ parameter space 
is constrained from all sides by indirect observations including structure formation, Ly Q forest, 
X-rays and phase space analysis, see section [3j 



2.5 Parameterization 

Adding k flavors of right handed neutrinos to the SM with three active neutrinos extends the 
parameter space of the model by 7k — 3 parameters. In the z/MSM k = 3, thus there are 18 
parameters in addition to those of the SM. These can be chosen as the masses rrii and Mj of 
the three active and three sterile neutrinos, respectively, and three mixing angles as well as three 
phases in each of the mixing matrices U v and Un that diagonalize m v and Mat, respectively. 

In the following we consider an effective theory with only two right handed neutrinos, which is 
appropriate to describe the generation of lepton asymmetries in scenarios I and II. After dropping 
N\ from the Lagrangian ([2]), the effective Lagrangian contains 11 new parameters in addition to 
the SM. 7 of them are related to the active neutrinos. In the standard parametrization they 
are two masses rrii (one active neutrino has a negligible mass), three mixing angles 0^, a Dirac 
phase 5 and a Majorana phase <fi. They can at least in principle be measured in active neutrino 
experiments. The remaining four are related to sterile neutrino properties. In the common Casas- 
Ibarra parametrization 56J] two of them are chosen as M 2 , M3. The last two are the real and 
imaginary part of a complex angle tu 0. The Yukawa coupling is written as 



F = U u yJmT g n^M. 



(15) 



where rn v 
given by 



diag 



diag(mi, 7ri2, 7713). For normal hierarchy of active neutrino masses (mi ~ 0) 1Z is 



/ 

1Z = cos(cu) sin(cu) 

\ — £sin(tu) £cos(tu) 

while for inverted hierarchy (7713 ~ 0) it reads 

cos(uj) sin(tu) 

1Z — I — £sin(cu) ^cos(oj) 





normal hierarchy 



(16) 



inverted hierarchy, 



(17) 



where £ = ±1. The matrix U u can be parameterized as 

U u = V 23 U s V 13 U_ s V l2 dmg(e ia ^ 2 ,e ia2 / 2 , 1) 



(18) 



'Note that F as a polynomial in z — e lco only contains terms of the powers z and 1/z. 
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sin 2 9i2 


sin 2 9i 3 


sin 2 9 23 


7.58 • 


2.35 • KT 3 


0.306 


0.021 


0.42 



Table 1: Neutrino masses and mixings as found in [571 /. We parameterize the masses m; according 
to mi = 0, m\ — m^i, m| = m^ tul + mf ol /2 for normal hierarchy and ml = mj tm — m^ ol /2, 
m| = JTiatm + m,3 = for inverted hierarchy. Using the values for 613 found more recently in 

pol . lE&j has no visible effect on our results. 





with U± s = diag(e Ti<5 / 2 , 1, e ±i5 ' 2 ) and 

Cl3 S13 

10 

-S13 C13 

where Cjj and Sjj stand for cos(9jj) and sin(0ij), respectively, and ai, 02 and 5 are the CP- 
violating phases. For normal hierarchy the Yukawa matrix F only depends on the phases ai 
and 5, for the inverted hierarchy, it depends on 5 and the difference a\ — 02- This is because 
iVi has no measurable effect on neutrino masses due to M\ <C M^^- In practice we will use the 
following parameters: two active neutrino masses rrii, five parameters in the active mixing matrix 
(three angles, one Dirac phase, one Majorana phase), the average physical sterile neutrino mass 
M = (M 1 + M 2 )/2 ~ M, the mass splitting AM. 

The masses and mixing angles of active neutrinos have been measured (the absolute mass 
scale is fixed as the lightest active neutrino is almost massless in scenarios I and II). We use the 



experimental values obtained from the global fit published in reference [53] in all calculations, 
which are summarized in table [TJ Shortly after we finished our numerical studies, the mixing 
angle 613 was measured by the Daya Bay [5H] and RENO [58| collaborations. The values found 
there slightly differ from the one given in [53], see also We checked that the effect of using 
one or the other value on the generated asymmetries in negligible, which justifies to use the self- 
consistent set of parameters given in table [TJ The remaining parameters can be constrained in 
decays of sterile neutrinos in the laboratory. 

It is one of the main goals of this article to impose bounds on them to provide a guideline for 
experimental searches. In order to identify the interesting regions in parameter space we proceed 
as follows. We neglect AM in (|15p . but of course keep it in the effective Hamiltonian introduced 
in section [U This is allowed in the region AM <C M, which we consider in this work. Unless 
stated differently, we always allow the CP-violating Majorana and Dirac phases to vary. We 
then numerically determine the values that maximize the asymmetry and fix them to those. In 
section [U where we study the condition i) for baryogenesis, we apply the same procedure to cu. 
On the other hand, the requirement iiV, necessary to explain SIdm in scenario I, almost fixes the 
parameter Recu to a multiple of 7r/2cj. In section [6] we therefore fix Recu = ir/2. 

The remaining parameters £, M, AM and Imcu contain a redundancy. For AM <C M 
changing simultaneously the signs of £, AM and Imtu along with the transformation Recu o 
7T — Recu corresponds to swapping the names of N2 and N3. To be definite, we always chose £ = 1 
and consider both signs of Imcu. Our main results consist of bounds on the parameters M, Imcu 
and AM. 



This is explained in section [276 
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For experimental searches the most relevant properties of the sterile neutrinos are the mass 
M ~ M and their mixing with active neutrinos. We therefore also present our results in terms 
of M, the physical mass splitting 5M and 



where G and U 2 are given by (|7|) and ([8]), respectively. U 2 measures the mixing between active 
and sterile species. SM and U 2 can, however, not be mapped on parameters in the Lagrangian 
in a unique way; there exists more than one choice of tu leading to the same U 2 . 

2.6 "Fine Tunings" and the Constrained vMSM 

In most models that incorporate the seesaw mechanism the eigenvalues of Mm are much larger 
than the scale of electroweak symmetry breaking. It is a defining feature of the i/MSM that all 
experimental data can be explained without introduction of such a new scale. In order to keep the 
sterile neutrino masses below the electroweak scale and the active neutrino masses in agreement 
with experimental constraints, the Yukawa couplings F have to be very small. As a consequence 
of this, the thermal production rates for lepton asymmetries are also very small unless they are 
resonantly amplified. In scenarios I and II this requires a small mass splitting between M2 and 
M3. This can either be viewed as "fine tuning" or be related to a new symmetry BE}. In the 
following we focus on these two scenarios, I and II. We do not discuss the origin of the small mass 
splitting here, but only list the implications EE 

Fermionic dispersion relations in a medium can have a complicated momentum dependence. In 
the following we make the simplifying assumption that all neutrinos have hard spacial momenta 
p ~ T and parameterize the effect of the medium by a temperature dependent quasiparticle 
mass matrix M/\r(T)f^l which we define as M^{T) 2 = H 2 — p 2 at |p| = p ~ T. Here H is 
the dispersive part of the temperature dependent effective Hamiltonian given in the appendix, 
cf. (|80p . The general structure of M/v(T) is rather complicated, but we are only interested in 
the regimes T < M (DM production) and T > Tew (baryogenesis) . Analogue to the vacuum 
notation in (fT0j) - (Tr2j) . we refer to the temperature dependent eigenvalues of Mn(T) as Mi(T) 
and M2(T), their average as M(T) and their splitting as 5M(T). Though Nj are the fields whose 
excitations correspond to mass eigenstates in the microscopic theory, the mass matrix Mjy(T) in 
the effective quasiparticle description is not necessarily diagonal in the iV^-basis for T / 0. The 
effective physical mass splitting 5M(T) depends on T in a non-trivial way. This dependence is 
essential in the regime M(T) 3> 5M(T), which we are mainly interested in. In principle also 
M(T) depends on temperature, but this dependence is practically irrelevant and replacing M(T) 
by M at all temperatures of consideration does not cause a significant error. 

There are three contributions to the temperature dependent physical mass splitting: The 
splitting AM that appears in the Lagrangian, the Dirac mass m,D{T) = Fv(T) that is generated 

12 As far as this work is concerned the sterile neutrino mass spectrum in the j^MSM follows from the requirement 
to simultaneously explain Qb and Qdm- It is in accord with the principle of minimality and the idea to explain 
new physics without introduction of a new scale (above the electroweak scale). In this work we do not discuss a 
possible origin of the mass spectrum and flavor structure in the SM and i/MSM, which to date is purely speculative. 
Some ideas on the origin of a low seesaw scale can be found in [6014711 ] . see also [72]. Some speculations on the 
small mass splitting have been made in 0, [jx| [4^ . Frll ] . A similar spectrum has been considered in a supersymmetric 




(20) 








for a discussion of the quasiparticle description. 
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by the coupling to the Higgs condensate and thermal masses due to forward scattering in the 
plasma, including Higgs particle exchangj^l. The interplay between the different contributions 
leads to non-trivial effects as the temperature changes. 



2.6.1 Baryogenesis 

For successful baryogenesis it is necessary to produce a lepton asymmetry of \i a ~ 10~ 10 at 
T > T + that survives until Tew an d is partly converted into a baryon asymmetry by sphalerons, 
see condition i). In this work we focus on scenarios I and II, in which only two sterile neutrinos 
^2,3 are involved in baryogenesis. In these scenarios baryogenesis is only possible if the physical 
mass splitting is sufficiently small (5M(T) <C M) and leads to a resonant amplification . On the 
other hand it should be large enough for the sterile neutrinos to perform at least one oscillation. 
Thus, baryogenesis is most efficient if it is of the same order of magnitude as the relaxation rate 
(or thermal damping rate) at T > T+, 

5M(T) ~ (T N )(T)u. (21) 

Here IV is the temperature dependent dissipative part of the effective Hamiltonian that appears 
in the kinetic equations given in section HI it is defined in appendix lA.3.2l and calculated in section 
14.21 It is essentially given by the sterile neutrino thermal width. However, (|2ip only provides 
a rule of thumb to identify the region where baryogenesis is most efficient. Numerical studies 
in section [5] show that the observed BAU can be explained even far away from this point, for 
M S> 5M(T) 3> Tjv(T). Thus, the mass degeneracy SM(T) <C M is the only serious tuning 
required in scenario II. In [3(3] it has been found that no such mass degeneracy is required in 
scenario III. 



2.6.2 Dark Matter Production 

In scenario I N\ dark matter has to be produced thermally from the primordial plasma [50(. In 
absence of chemical potentials, the resulting spectrum of N\ momenta has been determined in 



24i | . State of the art X-ray observations, structure formations and Ly a forest observations suggest 
that this production mechanism is not sufficient to explain Qdm because the required Aq mass 
and mixing are astrophysically excluded 21,[7(J. However, in the presence of a lepton chemical 



potential, the dispersion relation for active neutrinos is modified due to the MSW effect. If the 
chemical potential is large enough, this can lead to a level crossing between active and sterile 



neutrinos, resulting in a resonant amplification of the Aq production rate 52j. The full dark 
matter spectrum is a superposition of a smooth distribution from the non-resonant production 
and a non-thermal spectrum with distinct peaks at low momenta from resonant mechanism. In 
order to explain all observed dark matter by Aq neutrinos, lepton asymmetries |/x a | ~ 8 • 10~ 6 
are required at Tdm ~ 100 MeV [ljj. This is the origin of the condition ii) already formulated 
in section 12.41 Again the resonance condition (]2ip indicates the region where the asymmetry 
production is most efficient. For 1^, T_ <C Tew it imposes a much stronger constraint on the 
mass splitting than during baryogenesis because the thermal rates are much smaller. 

The asymmetries [i a can be created in two different ways, either during the freezeout of A2 3 
around T ~ T_ or in their decay at T ~ T^. During these processes we can use the vacuum 
value for v. As discussed in appendix I A. 3.11 the temperature dependence of 5M(T) is weak for 



We ignore the running of the mass parameters, which has been studied in [44 
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T < T_. The rates, on the other hand, still depend rather strongly on temperature, thus it 
is usually not possible to fulfill the requirement (I2TT) at T = T_ and T = simultaneously. 
Therefore one can distinguish two scenarios: the asymmetry generation is efficient either during 
freezeout (freezeout scenario) or during decay (decay scenario). On the other hand, (|21j) can be 
fulfilled simultaneously at T = T + and T = or at T = T + and T = T_ because at T = T + 
also the mass splitting depends on temperature. The strongest "fine tuning" requirement in the 
z^MSM is therefortEl 



5M(T + ) ~ (T N )u(T+) and 5M(T_) ~ (T N )u(T_) 
or 5M(T + ) ~ (rjv)/j(T+) and 5M(T d ) ~ (r^j^). 1 j 

From ([12]) it is clear that during the decay <5M(T d ) « 5M(T = 0) and 

5M " 2F( Re ( m ^ mD ) 3 3" Re ( m ^ mD ) 22 ) +AM (23) 
|5M| > l^{m) D m D ) 23 . (24) 



Fulfilling the resonance condition (121 p at low temperature requires a precise cancellation of the 
parameters in (j23|) and (|24p . both of which have to be fulfilled individually. The condition (|24p 
imposes a strong constraint on the active neutrino mass matrix ([5]). It can be fulfilled when real 
part of the off-diagonal elements is small. Note that this due to Q implies that Un is close to 
unity. This is certainly the case when the real part of complex angle cu in 1Z is a multiple of 
7r/2. In sections [6] and [7] we will focus on this region and always choose Recu = 7r/2. It should be 
clear that this is a conservative approach, since the production of lepton asymmetries can also be 
efficient away from the maximally resonant regions defined by (I22p . The lower bound ()23p can 
always be made consistent with (|2ip by adjusting the otherwise unconstrained parameter AM. 
At tree level this parameter is effectively fixed by 

AM = ~m ( Re H mD ) 33 " Re <H mD )J ± 5M ' (25) 

where the dependence of the RHS on AM is weak. The range of values for AM dictated by this 
condition is extremely narrow; it requires a tuning of order ~ 1CP 11 (in units of M). Quantum 



corrections are of order ~ rxii [4J], i.e. much bigger than 5M(T-). The high degree of tuning, 
necessary to explain the observed SIdm, is not understood theoretically. Some speculations can 



i explain tr 
be found in 

@, El H 123J]. However, the origin of this fine-tuning plays no role for the present 
work. 

In the following we will refer to the z^MSM with the condition Recu = ^ and the fixing of 
AM as constrained uMSM. Since the first term in the square root in (|12|) also depends on Recu, 
fixing this parameter exactly to a multiple of tt/2 usually does not exactly give the minimal 6M. 
However, it considerable simplifies the analysis, and deviations from such a value can in any case 
only be small due to the above considerations. 



It is in fact sufficient for baryogensis if 8M(T) ~ (Tn)ij(T) at some temperature T > T+ as long as some 
flavor asymmetries survive until Tew- The washout of the /i a typically becomes efficient around T+, but chemical 
equilibration can take long if one active flavour couples to the sterile neutrinos much weaker than the others. 
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3 Experimental Searches and Astrophysical Bounds 



The experimental, astrophysical and BBN bounds presented in this section and in the figures in 
sections [5][7] are derived under the premise that the mass and mixing of JV"i qualify it as a DM 
candidate, while iV^ are responsible for baryogenesis (scenarios I and II). Some of them loosen if 
one drops the DM requirement and considers the z/MSM as a theory of baryogenesis and neutrino 
oscillations only in scenario III. 



3.1 Existing Bounds 



A detailed discussion of the existing experimental and observational bounds on the z/MSM can 



be found in 2l| . Some updates that incorporate the effect of recent measurements of the active 
neutrino mixing matrix U u , in particular 613 7^ 0, have been published in [26l . 1271 . l29l ]. In the 



following we re-analyze all relevant constraints on the seesaw partners N23 from direct search 
experiments and BBN in the face of these experimental results. We also briefly review existing 
constraints on the dark matter candidate N\. 

As far as the known (active) neutrinos are concerned, the main prediction of the z^MSM 
is that one of them is (almost) massless. This fixes the absolute mass scale of the remaining 
two neutrinos [3]. Currently there is neither a clear prediction for the phases in U v in the 
z/MSM nor an experimental determination, though the experimental value for 813 55, 58| suggests 
that a measurement in principle might be possible. Regarding the sterile neutrinos, one has to 
distinguish between ^2,3 and N\. 



3.1.1 Seesaw Partners N2 and A3 

LHC The small values of Mj <C v in principle make it possible to produce them in the labo- 
ratory. However, the smallness of the Yukawa couplings F implies that the branching ratios are 
very small. Therefore the number of collisions (rather than the required collision energy) is the 
main obstacle in direct searches for the sterile neutrinos. In particular, they cannot be seen in 
high energy experiments such as ATLAS or CMS. It is therefore a prediction of the i^MSM that 
they see nothing but the Higgs boson. Vice versa, the lack of findings of new physics beyond the 
SM at the LHC can be viewed as indirect support for the model (though this prediction is of 
course relaxed if nature happens to be described by the i/MSM plus something else). 



Direct Searches The sterile neutrinos participate in all processes that involve active neutrinos, 
but with a probability that is suppressed by the small mixings U^. The mixing of to the SM 
is large enough that they can be found experimentally A number of experiments that allow 



to constrain the sterile neutrino properties has been carried out in the past [8C 



82j, |83|, NuTeV [8J, CHARM [85J], NOMAD |86j] and BEBC 



811 ] , in particular 
83] (see [si] for a 



CERN PS191 

review). These can be grouped into beam dump experiments and peak searches. 

Peak search experiments look for the decay of charged mesons into charged leptons (e 1 * 1 or 
fi^) and neutrinos. Due to the mixing of the active neutrino flavor eigenstates with the sterile 
neutrinos, the final state in a fraction of decays suppressed by (or is e^ + Nj (or fj^ + Ni). 
The kinematics of the two body decay can be reconstructed from the measured charged lepton, 
but the sterile flavor cannot be determined because of the Nj mass degeneracy. Therefore these 
experiments are only sensitive to the inclusive mixing U% defined in ([8]), where a is the flavor of 
the charged lepton. 
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In beam dump experiments, sterile neutrinos are also created in the decay of mesons, which 
have been produced by sending a proton beam onto a fixed target. A second detector is placed 
near the beamline to detect the decay of the sterile neutrinos into charged particles. Also in 
beam dump experiments, the sterile flavors cannot be distinguished. In this case, the expected 
signal is of the order U 2 U^ because creation and decay of the Nj each involve one active-sterile 
mixing. For instance, the CERN PS191 experiment constrains the combinations (U 2 ) 2 , (U 2 ) 2 , 



U 2 U 2 ance 



where 



E^M). ( 26 ) 



1 + 4 sin 2 6 W + 8 sin 2 W 1 - 4 sin 2 6 W + 8 sin 2 9 W 

, c^ = c T = . (27) 



This set differs from the quantities considered by the experimental group [82, |83j . It has been 



pointed out in [27j that the original interpretation of the PS191 (and also CHARM) data cannot 
be directly applied to the seesaw Lagrangian (pQ). The authors of [27j translate the bounds on 
active-sterile neutrino mixing published by the PS191 and CHARM collaborations into bounds 
that apply to the z/MSM and kindly provided us with their data. We use these bounds, along 
with the NuTeV constraints, as an input to constrain the region in the i/MSM parameter space 
that is compatible with experiments. 

Our results are displayed as green lines of different shade in the summary plots in figures 
[7j [13] and [14] in sections [5] and The different lines have to be interpreted as follows. Each 
shade of green corresponds to one experiment. For each experiment, there is a solid and a dashed 
line. The solid line is an exclusion bound. That means that there exists no choice of z/MSM 
parameters that leads to a combination of U 2 and M above this line and is consistent with table 
Q] and the experiment in question. In order to obtain the exclusion bound from an experiment 
for a particular choice of M we varied the CP- violating phases and Imcu0 We checked for each 
choice whether the resulting U 2 are compatible with the experiment in question. The exclusion 
bound in the M — U 2 plane is obtained from the set of parameters that leads to the maximal 
U 2 for given M amongst all choices that are in accord with experiment. The exclusion plots are 
independent of the other lines in the summary plots. The dashed lines (in the same shade as the 
exclusion plots) represent the bounds imposed by each experiment if the CP-violating phases are 
self-consistently fixed to the values that we used to produce the red and blue lines in the summary 
plots, which encircle the regions in which enough asymmetry is created to explain the BAU and 
DM. The NuTeV experiment puts bounds only on the mixing angle U 2 . This induces a much 
weaker constraint in the M — U 2 plane for inverted mass hierarchy than the other experiments. 
Our results differ from those of [22|. In the latter, the experimental constraints on the individual 



U 2 were directly reported in the M — U plane plotted in figure 3 of 22|. Moreover, only the 



PS191 exclusion bound was computed by distinguishing between mass hierarchies. 



Active Neutrino Oscillation Experiments The region below the "seesaw" line in figures [7] 
[T3l [TT1 and [14] is excluded because for the experimental values listed in table HJ there exists no 
choice of z^MSM parameters that would lead to this combination of M and U 2 . 

16 There are also constraints on U% which are, however, too weak to be of practical relevance. 
17 The mixings do not depend on Retu and the dependence on AM is negligible. 
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Big Bang Nucleosynthesis It is a necessary requirement that A^ have decayed sufficiently 
long before BBN that their decay products do not affect the predicted abundances of light ele- 
ments, which are in good agreement with observation. The total increase of entropy due to the 
A^2,3 decay is small, but the decay products have energies in the GeV range and even a small 
number of them can dissociate enough nuclei to modify the light element abundances. Since 
the sterile neutrinos are created as flavor eigenstates, they oscillate rapidly around the time of 
BBN. On average, they spend roughly half the time in each flavor state, and not the individual 
lifetimes of each flavor determine the relaxation time, but their average. This allows to estimate 
the inverse JV 2i 3 lifetime r by as r _1 ~ ^taT/v at T = 1 MeV. For r < 0.1s the decay products and 
all secondary particles have lost their excess energy to the plasma in collisions and reached equi- 
librium by the time of BBN0. We impose the condition r < 0.1s and vary all free parameters to 
identify the region in the M-U 2 plane consistent with this condition. The BBN exclusion bounds 
in figures [TJ [13] and HH represent the region in which no choice of z/MSM parameters exists that 
is consistent with table Q] and the above condition. Note that r" 1 ~ ^trTjv and the condition 
r < 0.1s are both rough estimates; the BBN bound we plot may change by a factor of order one 
when a detailed computation is performed. 



3.1.2 Dark Matter Candidate Ni 

The coupling of the DM candidate N\ is too weak to be constrained by any past laboratory 
experiment. However, different indirect methods have been used to identify the allowed region in 
the O^i — Mi plane. The possibility of sterile neutrino DM has been studied by many authors in 



the past, see |2ll. |89I| for reviews. In the following we summarize the most important constraints. 



As a decaying dark matter candidate N\ particles produce a distinct X-ray line in the sky 



that can be searched for. These pose an upper bound Mi < 3 — 4 keV (see e.g. |90h92|). If this 
were the only mechanism, the tension between these bounds would rule out the z/MSM as the 
common source of BAU, DM and neutrino oscillations. 

There are two different mechanisms for DM production in the i^MSM. The first one, common 
thermal (non-resonant) production leads to a smooth distribution of momenta. The second 
one, which relies on a resonance produced by a level crossing in active and sterile neutrino 
dispersion relations (see below), creates a highly non-thermal spectrum 1{J 5^1 . Observations of 



the matter distribution in the universe constrain the DM free streaming length. Without resonant 
production, the distribution reconstructed from Ly a forest observations suggests a lower bound 
on the mass Mi > 8 keV see also In combination with the X-ray bound, this 



would make resonant production necessary. In a realistic scenario involving both production 
mechanisms (\fi a \ ^ 10~ 5 ) this bound relaxes and has been estimated as Mi > 2 keV [79J. In our 
analysis we take these results for granted, though some uncertainties remain to be clarified, see 
section I3.2.H 

The DM production rate can be resonantly amplified by the presence of a lepton chemical 
potential in the plasma The resonance occurs due to a level crossing between active and 

sterile neutrino dispersion relations, caused by the MSW effect. This mechanism enhances the 
production rate for particular momenta as they pass through the resonance, resulting in a non- 



18 A more precise analysis of this condition for M < 140 MeV has been performed in [2£ 

19 There it has been assumed that baryonic feedback on matter distribution, see e.g. [96f | for a discussion, has 
negligible effect. 
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Figure 2: Different constraints on Ni mass and mixing. The blue region is e xclud ed by X-ray 
observations, the dark gray region M± < 1 keV by the Tremaine-Gunn bound ^9$t\l 0d /. The points 
on the upper solid black line correspond to observed Udm produced in scenario f in the absence of 
lepton asymmetries (for fj, a — 0) \19ij : points on the lower solid black line give the correct Qdm for 
|/ia| — 1-24 ■ 10 -4 , the maximal asymmetry we found. The region between these lines is accessible 
for < \/j, a | < 1.24 ■ 10~ 4 . We do not display bounds derived from Ly a forest observations because it 
depends on fi a in a complicated way and the calculation currently includes considerable uncertainties 



thermal DM momentum distribution that is dominated by low momenta and thus "colder" r°l 
Effectively, this mechanism "converts" lepton asymmetries into DM abundance, as the asymme- 
tries are erased while DM is produced. The full DM spectrum in the z^MSM is a superposition of 
the two components. The dependence on fj, a is, however, rather complicated. In particular, the 
naive expectation that the largest \[i a \, which maximized the efficiency of the resonant produc- 
tion mechanism, leads to the lowest average momentum ("coldest DM") is not true because ii a 
does not only affect the efficiency of the resonance, but also the momentum distribution of the 
produced particles. 

The N\ -abundance must correctly reproduce the observed DM density Qdm- This require- 
ment defines a line in the Mi — ^ Q |0 Q i| 2 -plane, the production curve. All combinations of Mi 
and Yl a \ ®ai\ 2 along the production curve lead to the observed DM abundance. Due to the res- 
onant contribution, the production curve depends on fi a . This dependence has been studied in 
fl9| , where it was assumed that jjL e = jj,^ = fM T . 



Finally, DM sterile neutrinos may have interesting effects for supernova explosions |10lMl05| |. 

Figure [2] summarizes a number of bounds on the properties of N±. The two thick black lines 
are the production curves for fi a = and |/x Q | = 1.24 • 10 -4 , the maximal asymmetry we found 
at T = 100 MeV in our analysis, see figure [12j The allowed region lies between these lines; above 
the fi a = line, the non-resonant production alone would already overproduce DM, below the 
production curve for maximal asymmetry fail to produce the required asymm etry for all 
choices of parameters. The maximal asymmetry has been estimated as ~ 7-10" 4 in 0,0, which 
agrees with our estimate shown in figure [12] up to a factor ~ 5. The corresponding production 



20 The results quoted in [97| demonstrate that the resonantly produced sterile is "warm enough" to change the 
number of substructure of a Galaxy-size halo, but "cold enough" to be in agreement with Ly a bounds ■ 
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curve is shown as a dotted line in figure [2j Our result is smaller and imposes a stronger lower 
bound on the iVi-mixing, which makes it easier to find this particle (or exclude it as the only 
constituent of the observed Qdm) in X-ray observations. However, though our calculation is 
considerably more precise than the previous estimate, the exclusion bound displayed in figure [12] 
still suffers from uncertainties of order one due to the issues discussed in appendix IA.4I and the 
strong assumption^! \i e = [l^ = /j, T , made in 0, [l9| to find the dependence of the production 
curve on the asymmetry. In order to determine the precise exclusion bound, the dependence of 
the production curve on individual flavor asymmetries has to be determined. 

The shadowed (blue) region is ex cluded by the non-observation of an X-ray line from N± decay 
in DM dense regions [7j, 



A lower limit on the mass of DM sterile neutrino 



below M\ ~ 1 keV can be obtained when applying the Tremai ne- Gu nn bound on phase space 



densities [981 ] to the Milky way's dwarf spheroidal galaxies [99j, [lOOf] . The Ly Q method yields 



similar bounds (not displayed in figure [2]). 



3.2 Future Searches 

3.2.1 Dark Matter Candidate Ni 

Indirect Detection The DM candidate N\ can be searched for astrophysically, using high 
resolution X-ray spectrometers to look for the emission line from its decay in DM dense regions. 
For details and refe renc es see the proposal submitted to European Strategy Preparatory Group 
by Boyarsky et al., 



U4]. 



Structure Formation Model- indep endent constraints on N% can be derived from consideration 
of dynamics of dwarf galaxies 98- lod ]. The existing small scale structures in the universe, such as 
galaxy subhalos, provides another probe that is sensitive to N\ properties because such structures 
would be erased if the mean free path of DM particles is too long. It can be exploited by comparing 
numerical simulations of structure formation to the distribution of matter in the universe that is 
reconstructed from Ly Q forest observations. However, the momentum distribution of resonantly 
produced N\ particles can be complicated [l9l . 52]. leading to a complicated dependence of the 
allowed mixing angle on the Ni mass and lepton asymmetries /x Q in the plasma [79I 92], A 
reliable quantitative analysis would involve numerical simulations that use the non-thermal N\ 
momentum distribution predicted in sc enar io I as input. While for Cold Darm Matter extensive 



studies have been performed, see e.g. [1151 ]. simulations for other spectra have only been done 



for certain benchmark scenarios; a model of Warm Dark Matter has been studied in [97]. 



Direct Detection As the solar system passes through the interstellar medium, the DM par- 
ticles Ni can interact with atomic nuclei in the laboratory via the n -\ 0* 1 suppressed weak in- 
teraction. This in principle opens the possibility of direct detection [U6fl. Such detection is 
extremely challenging due to the small mixing angle and the backgrou nd f r om solar and stellar 
active neutrinos. It has, however, been argued that it may be possible 117 . 1181 ]. 



21 Indeed, we find that this assumption does not hold in most of the parameter space. The asymmetries in 
individual flavors can be very different and even have opposite signs. The reason is that asymmetries generated 



at T > M are mainly "flavored" 
flavors can carry asymmetries. 



asymmetries, i.e. the total lepton number violation is small, but the individual 
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BBN The primordial abundances of light elements are sensitive to the number of relativistic 
particle species in the primordial plasma during BBN because these affect the energy budget, 
which determines the expansion rate and temperature evolution. Any deviation from the SM 
prediction is usually parameterized in terms of the effective number of neutrino species N e ff. 
At temperatures around 2 MeV, most N\ particles are relativistic. However, the occupation 
numbers are far below their equilibrium value, and the effect of the N% on N e ff is very small. 
Given the error bars in current measurements 1191 - 121 ]. the z^MSM predicts a value for N e ff 



that is practically not distinguishable from N e ff = 3. 

In principle the late time asymmetry in active neutrinos predicted by the i/MSM also affects 
BBN because the chemical potential modifies the momentum distribution of neutrinos in the 
plasma. However, the predicted asymmetry is several orders of magnitude smaller than existing 



bounds 4l[ and it is extremely unlikely that this effect can be observed in the foreseeable future. 



3.2.2 Seesaw Partners N 2 ,3 

The singlet fermions participate in all the reactions the ordinary neutrinos do with a probability 
suppressed roughly by a factor U 2 . However, due to their masses, the kinematic changes when 



an ordinary neutrino is replaced by Nj. The iv^ particles can be found in the laboratory |16[] 
using the strategies outlined in section [3.1. 1| which have been applied in past searches. 

One strategy, used in peak searches, is the study of kinematics of rare K, D, and B meson 
decays can constrain the strength of the Nj masses and mixings. This includes two body decays 
(e.g. — > /-i ± A r , K — > e ± A r ) or three-body decays (e.g. Kl,s — )• vr 1 * 1 + e T + ^2,3). The precise 
study of the kinematics is possible in O (like KLOE), charm, and B factories, or in experiments 
with kaons where the initial 4- moment umjs well known. For 3MeV < Mj < 414 MeV this 



possibility has recently been discussed in jl22l | 



The second strategy aims at observing the decay of the Nj themselves ( "nothing" — > leptons 
and hadrons) in proton beam dump experiments. The Nj are created in the decay K, D or B 
mesons emitted by a fixed target, into which the proton beam is dumped. The detector must 
be placed in some distance along the beamline. Several existing or planned neutrino facilities 
(related, e.g., to CERN SPS, MiniBooNE, MINOS or J-PARC), could be complemented by a 
dedicated near detector for these searches. Finally, these two stra tegies can be unified, so that 
the production and the decay occurs inside the same detector |l23 ]. 



For the mass interval M < mx, both strategies can be used. An upgrade of the NA62 
experiment at CERN would allow to search in the mass region below the Kaon mass mx- For 
m,K < M < mo it is unlikely that a peak search for missing energy at beauty, charm, and r 
factories will gain the necessary statistics. Thus, in this region the search for decays is 
the most promising strategy. Dedicated experiments using the SPS proton beam at CERN can 
completely explore the very interesting parameter range for M < 2 GeV. T his h as been outlined 



in detail in the European Strategy Preparatory Group by Gorbunov et al., (124 ]. 

An upgrade of the LHCb experiment could allow to combine both strategies. This would allow 
to constrain the cosmologically interesting region in the M — U 2 plane.. With existing or planned 
proton beams and B-factories the mass region between the D-mass and S-meson thresholds is 
in principle accessible, but such experiments would be extremely challenging. A search in the 
cosmologically interesting parameter space would require an increase in the present intensity of 
the CERN SPS bea m by two orders of magnitude or to produce and study the kinematics of more 



than 10 1U B -mesons [124J]. 
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4 Kinetic Equations 



Production, freezeout and decay of the sterile neutrinos are nonequilibrium processes in th e hot 



primordial plasma. We describe these by effective kinetic equations of the type used in [1251 ] 
and further elaborated in @, @, S, S3- These equations are similar to those commonly used to 
describe the propagation of active neutrinos in a medium. They rely on a number of assumptions 
and may require corrections when memory effects or off-shell contributions are relevant. These 
assumptions are discussed in appendix [A] We postpone a more refined study to the time when 
such precision is required from the experimental side. In the following we briefly sketch the 
derivation of the kinetic equations we use. More details are given in appendix lAl 



4.1 Short derivation of the Kinetic Equations 

We describe the early universe as a thermodynamical ensemble. In quantum field theory, any 
such ensemble - may it be in equilibrium or not - can be described by a density matrix p. The 
expectation value of any operator A at any time can be computed as 

(A)=tr(pA). (28) 

As there are infinitely many states in which the world can be, infinitely many numbers are 
necessary to exactly characterize p. These can either be given by all matrix elements of p or, 
equivalently, by all n-point correlation functions for all quantum fields. Either way, any practically 
computable description requires truncation. 

The leptonic charges can be expressed in terms of field bilinears, thus it is sufficient to con- 
centrate on the two-point functions. Instead of bilinears in the field operators themselves we 
consider bilinears in the ladder operators aj, a\ for sterile and a a , alt for active neutrinos. In 
principle there is a large number of bilinears, but only few of them are relevant for our purpose. 

For each momentum mode of sterile neutrinos we consider two 2x2 matrices formed by 
products of ladder operators a\aj, one for positive and one for negative helicity^l. Since M^r 
is diagonal in the iV/-basis, a\ai can be interpreted as a number operator for physical sterile 
neutrinos while a\aj with I J correspond to coherences. Nj are Majorana fields, but we can 
define a notion of "particle" and "antiparticle" by their helicity states. In the limit T 3> M, i.e. 
for a negligibly small Majorana mass term, the total lepton number (sum over a and j) defined 
this way are conserved. All other bilinears in the ladder operators for sterile neutrinos are either of 
higher order in F or quickly oscillating and can be neglected. Practically we are not interested in 
the time evolution of individual modes, but only in the total asymmetries. We therefore describe 
the sterile neutrinos by momentum integrated abundance matrices p^ for "particles" and pfj for 
"antiparticles" . The precise definitions are given in appendix lA.il 

The active leptons are close to thermal equilibrium at all times of consideration. This is 
because kinetic equilibration is driven by fast gauge interactions, while the relaxation rates for 
the asymmetries are of second order in the small Yukawa couplings F. We thus describe the active 
sector by four numbers^), the temperature and three asymmetries (one for each flavor, integrated 

22 This description is similar to the one commonly used in neutrino physics and could also be formulated in terms 
of "polarization vectors". 

23 It has been found in [l26| that mixing amongst the different SM lepton doublets can occur due to their coupling 
to the right handed neutrinos and affect leptogenesis. However, in the vMSM the Yukawa couplings F are too tiny 
to lead to a sizable effect. 
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over momentum). More precisely, the asymmetry in the SM leptons of flavor a is given by the 
difference between lepton and antilepton abundance, which we denote by p a , see (fl4"l) n . We 
study the time evolution of each flavor separately and find that they can differ significantly from 
each other. Following the steps sketched in appendix EJ one can find the effective kinetic "rate 
equations" 



. df>N 
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Here X = M/T, p eq is the common equilibrium value of the matrices pm and p^, H is the 
dispersive part of the effective Hamiltonian for sterile neutrinos that is responsible for oscillations 
and rates Tn, and form the dissipative part of the effective Hamiltonian. 

It is convenient to describe the sterile sector by p + and p_ , the CP-even and CP-odd deviations 
from equilibrium, rather than p^ and p^, 

px-ff9 = p+ -e=.. (32) 
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PN-P 

d29])- ([311 read 
.dp 



P+ + E T 



' dX 
.dp- 

'~dX 

. dp a 

1 dX 



[ReH, p + 
[ReH, p_ 
—iVJ^pa + S^, 



-{ReT N ,p + } + S_, 



-{ReT N ,p_} + S- 



with 



s„ 



- l ix + 2^ p - 



+ -{ImTjV,P-} 



1 



-p a ImT 



Ni 



2i[Imiz~, p + ] + {Imr^, p+} + ipaReT^ 



itr 



Re(f2)p- 



2tr 



Im(f£)^ 



(33) 
(34) 
(35) 

(36) 
(37) 
(38) 



Equations (|33j) - (j38|) are the basis of our numerical studies. 



4.2 Computation of the Rates 

The rates appearing in (f33"j) - (f38j) can be expressed as 

T N = (FaiKjRiT, M) aa 

a 

+F* I F aJ R M (T,M) aa ), (39) 
(f t)u (f%)u = x (F aI F*jR(T,M) aa 

-F*jF aJ R M (T,M) aa ), (40) 
T a L = x ((FF% a (R(T,M) aa + R M (T,M) aa )), (41) 

24 Note that the [i a here are abundances, not chemical potentials, which could alternatively be used to characterize 
the asymmetries. The relations between different characterizations of the asymmetries are given in appendix [Cl 
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Figure 3: Number of effective relativistic degrees of freedom g* as a function of temperature \l2^ ]. 

with no sum over a in (j40j) and (|4ip . The flavor matrices R and Rm are defined in appendix IA.3| 
see and (HU. Here F = FU N and 

dX M dT 2T 2 ' y ' 

where T 2 /M is the Hubble rate and M = M P (45) 1 / 2 /(47r 3 .^) 1 / 2 with the effective number of 
relativistic degrees of freedom computed in |l27l | and shown in figure El 

The flavor matrices R(T,M) and Rm(T,M) are almost diagonal since off-diagonal elements 
of Pafiip) include active neutrino oscillations, which are at least suppressed by rrii/T. We will 
always neglect the off-diagonal elements. R(T,M) and Rm(T, M) contain contributions from 
deca ys an d scatterings. In finite temperature field theory these can be associated with different 
cuts 



1281 ] through the Nj self-energy shown in figured! The scatterings keep the Nj in thermal 



equilibrium for T > T_. At T ~ T_ they become inefficient and the sterile neutrinos freeze out. 
Due to their small coupling they are long-lived, but unstable and decay at a temperature T^. For 
Tfi <C T_ , decay and freezeout are two separate processes and can be treated independently. This 
is the case in the interesting part of the z^MSM parameter space. 

4.2.1 Baryogenesis 

For T > v the SM fields are light and Mm is negligible. We can therefore neglect Rm(T, M) 
as well as the flavor dependence of R(T,M). Dispersion relations as well as relaxation rates are 
dominated by contributions from scatterings between vr and Higgs and lepton particles. The 
corresponding rates can be extracted from cuts through the diagrams shown in figure St). They 
have been computed in 0, H3] 

R(T, M) = 0.02^1 2x2 , R M (T, M) = for T > T EW . (43) 
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Figure 4: Contributions to the Ni self energies. Diagram a) dom inate s for T < v, diagram b) 
for T > v. Tjv is obtained fr om t he discontinuity of the diagrams '125], which can be computed 
by cutting it in various ways \l2$j . The gray self energy blobs indicate that dressed lepton and 
Higgs propagators have to be used. Cuts through them reveal a large number of processes, which 
are summarized in 123 / and appendix A of Q/. Recently it has been pointed out that current 
estimates suffer from an error 0(1) due to infrared and collinear enhancem ents at high temperature. 
Systematic approa ches to i nclude these effects can be found in \l3(X \l3l\ ] for T > M (relevant for 
baryogenesis) and \l3j\l3^] for M > T (relevant for late time asymmetries). We ignore this effect in 
our current study as it is comparable to other uncertainties in the kinetic equations and would only 
slightly change the results for the relevant regions in the vMSM parameter space. 
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Figure 5: The functions J? (s) (T,Af) and R,[ S f \T,M) for M = 1.2 GeV (darkest curve), M = 1.6 
GeV, M = 2 GeV, M = 2.5 GeV, M = 3 GeV, M = 3.5 GeV, M = 4 GeV (lightest curve). We 
would like to thank Mikko Laine for providing us numerical data. 



4.2.2 Dark Matter Production 

For Td <C which is the case in the interesting part of the z^MSM parameter space, freezeout 
and decay happen in different temperature regimes. At temperatures T > T_ the processes that 
keep the plasma in equilibrium are scatterings mediated by the weak interaction. Furthermore, 
the lepton masses are in first approximation negligible^!. Thus, R(T,M) and Rm(T,M) are 
proportional to the unit matrix and can be described by two scalar functions R( S \T,M) and 
R\j'{T,M). Around T_ ~ M, R^/(T,M) gives a small correction, which we account for in our 
analysis. In practice they have to be computed numerically. They are displayed in figure \5\ 

To be specific, in the high temperature limit, when all lepton masses are negligible, (|39p 
simplifies to 

T N ~ tU% ((F^F)*R {s \T, M) + F^FR$(T, M)) U* N at T ~ T_. (44) 
25 For some parameter choices this assumption can be violated for the r mass, introducing a small error. 
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In the low temperature regime, where R ~ Rm, one finds 

T N ~ t ^ (i^F* j) 4 D) (M) at T ~ T d . (45) 



The indices and indicate that the dominant contribution to the rate comes from scatterings 
or decays, respectively. The functions R^ (T, M) and R^p (T, M) can be obtained from the 
discontinuity of the Nj self energies shown in figure 2] at finite temperature. At T < M the 
different SM lepton masses become relevant. 

When the sterile neutrinos decay around Td <C M, the flavors are distinguishable. The 
decaying particle is non-relativistic and the density of the surrounding plasma low. For T = 
0, R(0,M) = Rm(0,M) take the same values and their elements are given by R(0,M) aa = 
Rm(0, M) aa = R^(M), where R&(M) are functions that can be computed from the vacuum 
decay rates for sterile neutrinos. The Nj can decay into various different final states, see appendix 
iDl There are leptonic and semi-leptonic channels, depending on the temperature either with 
quarks (before hadronization) or mesons (after hadronization) in the final state. Let V^ I ^ a be 
the rate at which Nj decays into a final state ip a of flavor a (e.g. v a qq or v a LpLp). Then 

i^)(M) = 2£%^, (46) 

where the sum runs over all possible final states that have flavor a. The factor 2 is due to the 
equal probabilities for decay into particles and antiparticles at tree level. The simple form of 
([4"8"]) is a result of the fact that the Yukawa couplings can be factored out of the corresponding 
amplitudes and the kinematics of N2 and A3 is the same due to their degenerate mass. Most 



of the rates required for our study have been computed in 16[], the remaining ones are given in 
appendix [Dj For T ~ 7^ with <C M the sterile neutrinos are non-relativistic and one can 
estimate 

Ad 

R(T <C M, M) aa ~ _i?(f)(M), (47) 

R M (T<.M,M) aa ~ Ez£M.rV>)(M), (48) 

h/ + p hi 

with E = (M 2 + p 2 ) 1 / 2 . Here p ~ T is the average sterile neutrino momentum. 
We therefore do not need all elements of the matrices R(T, M) and Rm(T, M) at all temperatures, 
but only two functions i?( 5 )(T, M) and R^\t, M) for T > T_ and three other functions r!£\m) 
for T <Td M. In practice we can simply add these contributions at all temperatures, though 
in principle we do not known the scattering contribution outside the range plotted in figure [5] and 
the decay contribution is obtained from vacuum rates. This is justified because for T > T_ our 
expressions for the decay rates are incorrect, but R& (M) <C R^\T,M). On the other hand 
R { a\M) > R ( ^{T,M) in the regime T <C M, T_, which is not covered by the data shown in 
figured! For < T < T_ our expressions for both, decay and scattering rates, are incorrect, 
but they are both smaller than the rate of Hubble expansion and have negligible effect. 

5 Baryogenesis from Sterile Neutrino Oscillations 

The BAU in the z^MSM is produced during the thermal production of sterile neutrinos Nj. This 
is in contrast to most other (thermal) leptogenesis scenarios, where decays and inverse decays 
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play the central role. The violation of total fermion number by the Majorana mass term Mm is 
negligible at Tew ^ M, but asymmetries in the helicity states of the individual flavors can be 
created. The sum over these vanishes up to terms suppressed by M/Tew, but because sphaleron 
processes only act on left chiral fields, the generated BAU can be much bigger. In this sense 
baryogenesis in the z/MSM can be regarded as a version of "flavored leptogenesis" . 

In this section we explore the part of the z^MSM parameter space where a BAU consistent 
with (|13p . i.e. r] ~ 10~ 10 , can be generated. We assume two sterile neutrinos iV^ participate in 
baryogenesis, as required in scenarios I and II. This assumption is motivated by the premise that 
N\ should be a valuable DM candidate, with masses and mixing consistent with astrophysical 
bounds. These require its Yukawa interaction to be too small to be relevant for baryogenesis, 
see section [3l In this sense, we consider the i/MSM as a model of both, baryogenesis and DM 
production, but are not concerned with the DM production mechanism, which is discussed in 
the following section [6l This corresponds to scenario II. The requirement to explain Qdm only 
enters implicitly, as we demand the N\ mass and mixing to be consistent with astrophysical 
observations. If one completely drops the requirement to explain the observed DM and study 
the i^MSM as a theory of baryogenesis and neutrino oscillations only (scenario III) , the resulting 
bounds on the parameters weaken considerably. In particular, it was found in [30( that no mass 
degeneracy between the sterile neutrino masses is needed. 

We extend the analysis performed in [22I], but take into account two additional aspects. First, 
we use the non-zero value for the active neutrino mixing angle 813 given in tableQ] which brings in 
a new source of CP-violation through the phase 8. Second, we include the contribution from the 
temperature dependent Higgs expectation value v(T) to the effective Hamiltonian, coming from 
the real part of the diagram in figure Uk). It is relevant for temperatures close to the electroweak 
scale. 

We solve numerically the system of equations f|33j) -f)38 [) to find the lepton asymmetries at 
T ~ Tew 1 assuming that there is no initial asymmetry. The effective Hamiltonian is given by 
(|84p and (|39p - ([4ip with (|43p . We fix the active neutrino masses and mixing angles according 
to table Q] and choose the phases 5, a% and 02 as well as Recu to maximize the asymmetry. 
Interestingly, for normal hierarchy of neutrino masses, the value of Recu that maximizes the 
asymmetry is close to ^, as required in the constrained i/MSM. This allows to identify the region 
in the remaining three-dimensional parameter space consisting of M, AM and Imcu where an 
asymmetry > 10~ 10 can be created. Deep inside this region, the asymmetry generated for this 
choice of phases can be much too large, but it can always be reduced by choosing different phases. 
Thus, any choice of M, AM and Imcu inside this region can reproduce the observed BAU. 

In practice it is difficult to find the phases that maximize the asymmetry in each single 
point, as we are dealing with a seven-dimensional parameter space. However, the analysis can be 
simplified. First, the choice of phases that maximize the asymmetry practically does not depend 
on AM because the dependence of the Yukawa coupling (|15p on AM is very weak. Second, our 
numerical studies reveal in most of the parameter space Imcu is the main source of CP-violation. 
The other phases have comparably little effect on the final asymmetry, except for the region 
around Imcu = 0. Surprisingly, the values for 8, a\ and a<i that maximize it vary only very little 
and are always close to zero. One possible interpretation is that Imcu provides the main source for 
the asymmetry generation, while 8, a\ and (X2 contribute stronger to the washout. However, due 
to the various different time scales involved we cannot extract a single CP-violating parameter at 
this point, which is commonly used in thermal leptogenesis scenarios to study such connections 
analytically. The above seems to be valid everywhere except in the region Imcu ~ 0, where 8, a\ 
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and «2 are the only sources of CP-violation. 

We therefore split the parameter space into two regions. In the region 0.5 < e lmw < 1.5 we 
chose «2 = 7r, 5 = 0, Recu = j^tt for normal hierarchy and «2 — &i = tt, 5 = tt, Retu = |-7r 
for inverted hierarchy. Everywhere else we chose a?2 = -Tp- and 5 = 7r, Recu = ^7r for normal 
hierarchy and c*2 — ol\ = y^7r, 5 = ^7r, Retu = |-7r for inverted hierarchy. Note that for normal 
hierarchy F only depends on 02 and 6, while for inverted hierarchy it depends on 02 — ol\ and 6 
because one neutrino is massless. In order to determine v(T) one nee ds to fix t he Higgs mass ran ■ 
We used the value mu = 126 GeV suggested by recent LHC data |l35l . Il36l |. corresponding to 
electroweak scale of Tew ~ 140 GeV. This is consistent with the i/MSM being a valid description 
of nature up to the Planck scale 1371 ] . 



We present our results in figure El which shows the allowed region in the AM — Imcu plane for 
several masses M. The lines correspond to the exact observed asymmetry, inside more asymmetry 
is generated. As pointed out above, any point inside the lines is consistent with observation 
because the asymmetry can be reduced by choosing different phases. Figure [6] shows that even 
for small masses around 10 MeV enough asymmetry can be created. However, for small masses 
the CP-violation contained in 5, a\ and 02 is not sufficient, and the allowed region consists of two 
disjoint parts that are separated by the Imcu ~ region. The area of these increases with M. For 
masses of a few GeV, the CP-violation from 5, ol\ and 02 alone is sufficient and the regions join. 
Interestingly, there appear to be mass-independent diagonal lines in the AM — e Imcu plane that 
confine the region where enough asymmetry can be generated. We currently have not understood 
the origin of these lines parametrically. The inverted hierarchy generally allows to produce more 
asymmetry than the normal hierarchy. There is an approximate symmetry between regions with 
positive and negative Imcu. It would be exact when simultaneously changing £ and is related 
to the symmetry of the Lagrangian under exchange of N2 and N3. As expected, these results 
are close to those obtained in [22|, which provides a good consistency check. The slightly bigger 
asymmetry is due to the additional source of CP-violation for 813 7^ 0. 

For experimental searches, the most relevant parameters are the mass M and the mixing 
between active and sterile neutrinos. In figure [7] we translate our results into bounds on the 
flavor independent mixing parameter U 2 defined in equation (|20p . Using the results displayed 
in figure (H we chose 5M to maximize the asymmetry and find the region in the U 2 — M plane 
within which baryogenesis is possible. The plot has to be read as follows: For each point in the 
region between the blue lines there exists at least one choice of z^MSM parameters that allows for 
successful baryogenesis. The plots in figure [7] are similar to the ones of figure 3 in ^2^, but the 
allowed region is slightly bigger due to the effect the new source of CP- violation for 613 7^ 0. 

The constraints on the mixing angle U 2 shown in figure [7] can be translated into constraints 
on the neutrino lifetime r _1 ~ ^trr^ (at T = 1 MeV) shown in figure The plot in figure [8] is 
similar to the ones of figure 4 in 



6 Late Time Lepton Asymmetry and Dark Matter Production 

The lepton asymmetry at temperatures of a few hundred MeV is of crucial importance for the dark 
matter production in scenario I. Resonant dark matter production requires a lepton asymmetry 
\n a \ ~ 8- 10~ 6 in the plasma, much larger than the baryon asymmetry. The details of this process 
have been outlined in 0] . Here we are not concerned with the dark matter production itself, but 
with the mechanisms that generate the required lepton asymmetry. This asymmetry must come 
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Figure 6: Values of AM and Imtu that lead to the observed baryon asymmetry in scenarios I and II 
for different sterile neutrino masses M = 10, 100 MeV, 1 and 10 GeV. The blue (shortest dashed) 
line corresponds to M = 10 MeV, red (short dashed) to M = 100 MeV, brown (long dashed) to 
M = 1 GeV and green (longest dashed) to M = 10 GeV. The phases that maximize the asymmetry 
differ significantly for Imtu » and away from that region. In the region 0.5 < e Imco < 1.5 we chose 
ct-2 = it, 8 — 0, Recu = y^7r for normal hierarchy and 02 — o.\ — 7r, 8 = ir, Recu = |7r for inverted 
hierarchy. Everywhere else we chose 8 = ^jtt, Retu = \n for normal hierarchy and 02 — ai = y^vr, 
8 = go 7T, Recu = |7r for inverted hierarchy. The upper panel shows the results for normal hierarchy, 
the lower panel for inverted hierarchy. 

from a source that is different from that of the baryon asymmetry because reach chemical 
equilibrium at some temperature T+ < Tew and the asymmetry in the leptonic sector is washed 
out (while the baryon asymmetry remains as sphalerons are inefficient at T < Tvmv r^l. There 
are two distinct mechanisms that contribute to the late time asymmetry, the freezeout of at 

26 It has been suggested that some asymmetry may be preserved in magnetic fields down to temperatures T < T_ 
due to the chiral anomaly [49(] . Here we take the most conservative approach and do not take into account this 
possibility. 
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Figure 7: Constraints on the N2.3 masses A/2,3 — M and mixing U 2 = tr(6^6) from baryogenesis in 
scenarios I and II; upper panel - normal hierarchy, lower panel - inverted hierarchy. In the region 
between the solid blue "BAU" lines, the observed BAU can be generated. The regions below the solid 
black "seesaw" line and dashed black "BBN" line are excluded by neutrino oscillation experiments 
and BBN, respectively. The areas above the green lines of different shade are excluded by direct 
search experiments, as indicated in the plot. The solid lines are exclusion plots for all choices of 
vMSM parameters, for the dashed lines the phases were chosen to maximize the BAU, consistent 
with the blue lines. 



T ~ T_ and their decay at T ~ 

The requirement that these two mechanisms produce enough asymmetry put severe con- 
straints on the parameters of the model, described in section 12.61 The value of Recu is fixed to 
values near ir/2. The mass splitting AM is limited to a very narrow range by equation (|25|) . 
Therefore we will use the mass splitting in vacuum 5M instead of AM as a free parameter in 
the following. All experimentally known parameters are fixed to the values given in table [TJ The 
phases 5, oi\ and «2 are chosen to maximize the asymmetry. As in section [5] we observe that in 
most of the parameter space Imcu is the main source of CP-violation. We again find that it is 
convenient to split the parameter space into the region 0.5 < e lmw < 1.5 and the complement. 
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Figure 8: Constraints on the N2,3 masses Mi,z ~ M and lifetime r _1 ~ |trrjv|T=iMcV from baryo- 
genesis in scenarios I and II. In the region between the blue "BAU" lines, the observed BAU can be 
generated (solid line - normal hierarchy, dotted line - inverted hierarchy). The region above the solid 
black "BBN" line is excluded by BBN. 

For normal hierarchy we chose the phases a 2 = § > 8 = I 71 " in the region 0.5 < e lmw < 1.5 and 
q.2 = § and S = everywhere else. For inverted hierarchy we chose 02 — c^i = 5 and 5 = |tt in 
the region 0.5 < e lmw < 1.5 and «2 — Q-i = 0, 5 = yjj7r everywhere else. Note that for normal 
hierarchy F only depends on «2 and <5, while for inverted hierarchy it depends on 02 — ot\ and 5 
because one neutrino is massless. We then study the parameter space spanned by M, 5M and 
Imtu. 

As in section [H we use the kinetic equations (]33p -(|38 p in order to calculate the lepton asym- 
metries as a function of T. The effective Hamiltonian is calculated from (|87p and (|39p -(|4ip with 
(|44p - (l48p . We impose thermal equilibrium with vanishing chemical potentials as initial condition 
at a temperature T > T_ and look for the parameter region where ^ a \ ^a\ > 8 • 10~ 6 at T = 100 
MeVllI 

The results are shown in figures [9] and [TUJ The required asymmetry can be created when 
the sterile neutrinos have masses in the GeV range. For small masses of M = 2 — 4 GeV the 
CP violation contained in ai, 02 and 5 alone is not sufficient for normal hierarchy and barely 
sufficient for inverted hierarchy; a non-zero Imtu is required and the allowed region consists of two 
disjoint parts along the Imtu axis which are separated by the Imtu ~ region. For larger masses 
(M > 7 GeV for normal hierarchy, M > 4 GeV for inverted hierarchy ), the regions merge, but 
Imtu continues to be the most relevant source of CP violation in most of the parameter space. In 
addition, one can also observe disjoint regions along the 5M axis. These can be identified with 
the decay scenario and freezeout scenario. In the upper part of the figures, the asymmetry is 
mainly created during the freezeout of ^2,3, in the lower part during the decay. At T_, Tjy has 
considerably larger entries than at T^. Thus, the resonance condition (I2ip requires a smaller mass 
splitting in the decay scenario. For larger masses, both regions merge. However, freezeout and 
decay are always two separated processes, i.e. T_ S> T^. As in figure (U there is an approximate 
symmetry under a change of sign for Imtu, which is related to the symmetry of the Lagrangian 

27 We solve the kinetic equations down to T — 50 MeV in order to avoid numerical artifacts at the boundary. 
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under exchange of N2 and N% and becomes exact when also changing £. 



10-" 




0.1 0.2 0.5 1.0 2.0 5.0 10.0 20.0 




0.1 0.2 0.5 1.0 2.0 5.0 10.0 20.0 




0.1 0.2 0.5 1.0 2.0 5.0 10.0 20.0 




0.5 1.0 2.0 

JmM 



5.0 10.0 20.0 



Figure 9: VaJues of 5M and Imcu that lead to the lepton asymmetry required for dark matter 
production in scenario I for different singlet fermion masses, M — 2.5, 4, 7 and 10 GeV and for 
normal hierarchy. The upper left panel corresponds to M = 2.5 GeV, the upper right panel to 
M = 4 GeV, the lower left panel to M = 7 GeV and the lower right panel to M = 10 GeV. The 
phases that maximize the asymmetry differ signihcantly for Imiv ~ and away from that region. We 
chose the phases 0:2 = -f , 8 = ~7T in the region 0.5 < e ImuJ < 1.5 and ct2 = \ and 5 = everywhere 
else. 



For experimental searches for sterile neutrinos, the most relevant parameters are the mass M 
and the mixing between active and sterile species. As in section [5j we translate our results for 
the parameters in the Lagrangian into bounds on the mass and mixing. For each mass, we chose 
SM in a way that maximizes the allowed region in the U 2 — M-plane. The results are shown in 
figure [TTJ 

Finally, we estimate the maximal asymmetry that can be generated at T ~ 100 MeV as a 
function of M by its largest value within the data files we used to create figures and [TUJ The 
maximal asymmetry allows to impose a lower bound on the N\ mixing; bigger lepton asymmetries 
make the resonant DM production more efficient and allow for sm aller N\ mixing, displayed in 
figure [2j Furthermore, the maximal \fx a \ is of interest because in 1381 ] it was pointed out that 
a large lepton asymmetry may lead to a first order phase transition during hadronisation. The 
maximal asymmetries we found are shown in figure Q2J For both hierarchies they remains well 
below cosmological bounds (see 0]) at all masses of consideration and are about a factor 5 
smaller than the value 7 • 10 -4 estimated in 0]. However, given the uncertainties summarized in 
appendix I A.41 they can easily change by a factor 0[1}. 
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Figure 10: Values of 8M and Imiv that lead to the lepton asymmetry required for dark matter 
production in scenario I for different singlet fermion masses, M = 2.5, 4, 7 and 10 GeV and inverted 
hierarchy. The upper left panel corresponds to M = 2.5 GeV, the upper right panel to M = 4 GeV, 
the lower left panel to M = 7 GeV and the lower right panel to M — 10 GeV. The phases that 
maximize the asymmetry differ significantly for Imw ~ and away from that region. We chose 
Oi'z — on = t and 8 — f 7r in the region 0.5 < e Irnw < 1.5 and ct2 — ct\ = 0, 8 = everywhere else. 
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Figure 11: Constraints on the N2,?, masses Af2,3 — M and mixing U 2 = tr(9^6) in scenario f. The 
lepton asymmetry at T — 100 MeV can be large enough that the resonant enhancement of Ni 
production is sufficient to explain the observed Qdai inside the dashed blue and red lines for normal 
and inverted neutrino mass hierarchy, respectively. The regions below the "seesaw" lines are excluded 
by neutrino oscillation experiments for the indicated choice of hierarchy. 
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Figure 12: Estimate of the maximal lepton asymmetry that can be created in scenario I around 
T = 100 MeV. The blue small dashed line line corresponds to normal hierarchy, the red long dashed 
line to inverted hierarchy. This plot was generated using the maximal values found in the data files 
used to produce Gsures W\ and 1 1 01 

7 DM, BAU and Neutrino Oscillations in the z/MSM 

In the previous sections and [U] we have studied independently the conditions for successful 
baryogenesis on one hand and sufficient dark matter production on the other. The most interesting 
question is of course in which part of the z/MSM parameter space scenario I can be realized, 
i.e. both can be achieved simultaneously. This region cannot be found by simply superposing 
the figures from the previous sections because the phases that maximize the asymmetry are 
different for T > Tew an d T < T_. The requirement to produce enough DM imposes the 
stronger constraint. We therefore fix the CP-violating phases in a way that is consistent with 
X^a l/^ol > ^ • 10~ 6 at T = 100 GeV in some significant region in the M — C/ 2 -plane. We then 
check for which combination of M and U 2 the correct BAU is created by these phases. 

We start with the phases used in figure [TTl which were chosen to maximize the area in the 
M— C/ 2 -plane where |/x a | > 8T0~ 6 at T = 100 MeV. The result is shown in figure[13j The blue 
line corresponds to the points where the asymmetry at Tew corresponds to the observed BAU. 
While the requirement to produce enough DM only imposes a lower bound on the asymmetry at 
100 MeV, the value of the BAU is known to be given by (|13p . i.e. has a fixed value. Thus, only 
the points on the blue line that lie within the region encircled by the red line (DM region) form 
the allowed parameter space. The shape of the blue BAU line can be modified by changing the 
phases, see figure [U but this also changes the shape of the red line (DM region). Solving the 
kinetic equations for different phases reveals that the BAU line can be brought to most points 
within the DM region. This region therefore gives a good estimate of the allowed parameter 
space. 

The constraints derived on the mixing angle U 2 are translated into constraints on the neutrino 
lifetime r" 1 ~ ±trT N (at T = 1 MeV) in figure [EJ 

In the plots of figure [T¥| there are two regions where the 'BAU' and 'DM' lines are close, 
leading to the successful baryogenesis and dark matter production. One is near the seesaw line 
and the other is for higher mixing. These regions are easier to identify in the Imtu — M plane 
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Figure 13: Constraints on the ^2,3 masses M2,3 — M and mixing U 2 = tr(6' r 9) in the constrained 
vMSM (scenario I); upper panel - normal hierarchy, lower panel - inverted hierarchy. In the region 
between the solid blue "BAU" lines, the observed BAU can be generated. The lepton asymmetry 
at T — 100 MeV can be large enough that the resonant enhancement of Ni production is sufficient 
to explain the observed Qdm inside the solid red "DM" line. The CP-violating phases were chosen 
to maximize the asymmetry at T = 100 MeV. The regions below the solid black "seesaw" line and 
dashed black "BBN" line are excluded by neutrino oscillation experiments and BBN, respectively. 
The areas above the green lines of different shade are excluded by direct search experiments, as 
indicated in the plot. The solid lines are exclusion plots for all choices of vMSM parameters, for the 
dashed lines the phases were chosen to maximize the late time asymmetry, consistent with the red 
line. 




shown in figure [T6j The baryon asymmetry almost vanishes for Imcu really close to 0, but this 
is not the case for dark matter production. Therefore, there is a region near Imcu = which 
produce the right amount of baryon asymmetry and enough dark matter. This is the region 
where the blue 'BAU' line is inside the red 'DM' line in figure [TBI For large value of |Imcu|, there 
also are regions where the two constraint are close. 
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Figure 14: Same as figure \13\ but with a different set of CP-violating phases. The plot illustrates 
how the "BAU" line moves when the phases are changed; it can cross through points deep inside 
the maximal "BAU" region while the phase still allow for DM production. The sign of the BAU is 
opposite in the two disjoint "BAU" regions. 
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Figure 15: Constraints on the N2.3 masses Mi,z — M and lifetime t~ ~ |trFjv|T=iMeV in the 
constrained vMSM (scenario I). In the region between the blue "BAU" lines, the observed BAU 
can be generated. The lepton asymmetry at T = 100 MeV can be large enough that the resonant 
enhancement of Ni production is sufficient to explain the observed Qdm inside the red "DM" line. 
The CP-violating phases were chosen to maximize the asymmetry at T — 100 MeV. Solid lines - 
normal hierarchy, dotted lines - inverted hierarchy. 

8 Conclusions and Discussion 

We tested the hypothesis that three right handed neutrinos with masses below the electroweak 
scale can be the common origin of the observed dark matter, the baryon asymmetry of the universe 
and neutrino flavor oscillations. This possibility can be realized in the i/MSM, an extension of the 
SM that is based on the type-I seesaw mechanism with three right handed neutrinos Nj. Center- 
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Figure 16: Constraints on the A^,3 masses Mifi — M and parameter Imu> in the constrained vMSM 
(scenario I); upper panel - normal hierarchy, lower panel - inverted hierarchy. In the region between 
the solid blue "BAU" lines, the observed BAU can be generated. The lepton asymmetry at T = 100 
MeV can be large enough that the resonant enhancement of Ni production is sufficient to explain 
the observed CIdm inside the solid red "DM" line. The CP-violating phases were chosen to maximize 
the asymmetry at T — 100 MeV. 

piece of our analysis is the study of sterile and active neutrino abundances in the early universe, 
which allows to determine the range of sterile neutrino parameters in which DM, baryogenesis 
and all known data from active neutrino experiments can be explained simultaneously within the 
z^MSM. We combined our results with astrophysical constraints and re-analyzed bounds from 
past experiments in the face of recent data from neutrino oscillation experiments. We found that 
all these requirements can be fulfilled for a wide range of sterile neutrino masses and mixings, see 
figures [T3l [TH in section [7J In some part of this parameter space, all three new particles may be 
found in experiment or observation, using upgrades to existing facilities. 

This is the first complete quantitative study of the above scenario (scenario I), in which no 
physics beyond the z^MSM is required. We found that the i/MSM can explain all experimental 





37 



data if one sterile neutrino (Aq), which composes the dark matter, has a mass in the keV range, 
while the other two (A^) have quasi-degenerate masses in the GeV range. The heavier particles 
A^2,3 generate neutrino masses via the seesaw mechanism and create flavored lepton asymmetries 
from CP-violating oscillations in the early universe. These lepton asymmetries are crucial on two 
occasions in the early universe. On one hand they create the BAU via flavored leptogenesis. One 
the other hand they affect the rate of thermal DM production via the MSW effect. The second 
point allows to derive strong constraints on the A/2,3 properties from the requirement to explain 
the observed Qdm by thermal Aq production, see section El This can be achieved by resonant 
production, caused by the presence of lepton asymmetries in the primordial plasma at T ~ f 00 
MeV. The required asymmetries can be created when A^ are heavier than 1 — 2 GeV and the 
physical mass splitting between the A2 and A3 masses is comparable to the active neutrino mass 
differences. This can be achieved in a subspace of the z^MSM parameter space that is defined by 
fixing two of the unknown parameters (the Majorana mass splitting AM and a mixing angle Recu 
in the sterile sector). This choice, in which scenario I can be realized, is dubbed "constrained 
z/MSM". 

We also studied systematically how the parameter constraints relax if one allows Aq DM to 
be produced by some unspecified mechanism beyond the z^MSM (scenario II), see section [5j In 
this case the strongest constraints come from baryogenesis and the required mass degeneracy is 
much weaker, AM/M < 10~ 3 . We found that successful baryogenesis is possible for A9.3 masses 
as low as 10 MeV. These results are based on an extension of the analysis performed in [221 ] that 
accounts for a non-zero value of the neutrino mixing angle 813 and a temperature dependent Higgs 
expectation value. While the low mass region is severely constrained by BBN and experiments, 
the allowed parameter space becomes considerably bigger for masses in the GeV range. Detailed 
results for the allowed sterile neutrino masses and mixings are shown in figures [6] - [71 

If one completely drops the requirement that DM is composed of Aq and considers the i/MSM 
as a theory of baryogenesis and neutrino oscillations only (scenario III), no degeneracy in masses 
is required. Note that this also implies that no degeneracy is required in scenario II if more than 
three right handed neutrinos are added to the SM. 

For masses below 5 GeV, the heavier sterile neutrinos can be searched for in experiments 
using present day technology. This makes the z/MSM one of the few truly testable theories of 
baryogenesis. The parameter space for the DM candidate Aq is bound in all directions, see figure 
[2j and can be tested using observations of cosmic X-rays and the large scale structure of the 
universe. Since the model does not require new particle physics up to the Planck scale to be 
consistent with experiment, the hierarchy problem is absent in the z/MSM. 

We conclude that neutrino physics can explain all confirmed detections of physics beyond the 
standard model except accelerated cosmic expansion. 
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A Kinetic Equations 



In the following we sketch the derivation of the kinetic equations (|29p -(|3ip. Our basic assumptions 
can be summarized as follows: 

1) Coherent states containing more than one Nj quantum are not relevant. Their contributions 
are suppressed by additional powers of the small mixing 8 a j. Processes involving one 
sterile neutrino include decays of Nj particles, their scatterings with SM particles and 
flavor oscillations. 

2) Screened one-particle states are the only relevant propagating neutrino degrees of freedom. 
In particular, we do not consider any collective excitations, which are infrared effects and 
only give a small contribution when the typical neutrino momenta are hard ~ T. 

3) The interactions that keep the SM fields in equilibrium act much faster than interactions 
involving Nj at all times due to the smallness of F. 

4) T_ > Trf, i.e. the lifetime of the Nj is sufficiently long that freezeout and decay are two 
well-separated events. This is the case in the parameter space we study. 

5) The typical momentum of Nj particles is p ~ T even when they are out of equilibrium. 
This is justified because they are produced from a thermal bath and freeze out from a 
thermal state, hence their distribution functions should mimic those of kinetic equilibrium 
even when out of equilibrium. 

6) We neglect the effect of the on the time evolution of the entropy (or temperature). 
This is justified as their contribution to the total entropy and energy densities are always 
small. 

7) We neglect t he ef fect of the lepton asymmetry on hadronization. This aspect has e.g. been 



discussed in [1381 ] 



A.l How to characterize the Asymmetries 

The leptonic charges that we are interested in can be expressed in terms of field bilinears. This 
and 1) imply that this is sufficient as we only need to deal with a reduced density matrix, in which 
all states including more than one sterile neutrino have been removed by partial tracing. Instead 
of bilinears in the field operators themselves we consider expectation values of bilinears in the 
ladder operators ai, a\ for sterile and a a , a a for active neutrinos. To be explicit, we decompose 
Nj a@ 

Nl = I ^E^Kp^^^P'^ + ^p^H^P'*))- ( 49 ) 

Here p is momentum, s the spin index and u, v are the usual plane wave solutions to the Dirac 
equation, 

(y-MjK )P = (j/+M I )vl p = 0. (50) 

28 The sterile neutrinos may be described as Weyl, Dirac or Majorana spinors. Here we have chosen to write vr 
as a right chirality four-spinor and pulled the Pr out of the definition Q so that the Ni are Majorana spinors. 
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We will in the following always assume that the spacial momentum is directed along the z- 
axis, which is also the axis of angular momentum quantization. We chose the convention that 



hup = (— l) s+1 «p while hup = (— l) s f p , where h is the helicity matrix 



1 Pi i 5 

-—7 7 7 

2 p 



2|P| 1 



(51) 



All relevant matrix elements of the density matrix p can be identified with expectation values 
of bilinears in the ladder operators. Because of this the matrix p of bilinears, to be defined in 
(|52p . is often referred to as "density matrix" (rather than p itself). In principle there is a large 
number of such bilinears. A complete characterization of the system requires knowledge of all 
their expectation values at all times. However, it can be simplified dramatically, and for our 
purpose it will be sufficient to follow the time evolution of two 2x2 matrices pn and pjy and 
three chemical potentials. 

The only term in (JT]) that violates lepton number is Mm- For T ^ M, it is negligible and 
lepton number is approximately conserved. There is no total lepton asymmetry at Tew in the 
z/MSM, but there can be asymmetries of opposite sign for fermions with different chirality. Baryo- 
genesis occurs because sphalerons only couple to left handed fermions. As far as the (Majorana) 
neutrinos are concerned, the two helicity states act as "particle" and "antiparticle" . Terms con- 
taining two creation or two annihilation operators such, as (ajaj) or (a„a^), can be related to 
processes that violate lepton number and are suppressed at T > M. For T < M they in principle 
could contribute, but the leading order contribution in the Yukawa coupling F to the correspond- 
ing rates ^(ajaj) etc@ oscillate fast. We therefore only consider terms that contain exactly one 
creation and one annihilation operator. Since only two of the sterile neutrinos are relevant here, 
these form a 10 x 10 matrix that can be written as 

/ Pnn Pnn Pnl Pnl \ 

Pnn Pnn Pnl Pnl 

Pln Pln Pll Pll 

V Pln Pln Pll PLL ) 



(52) 



with 



(pnn)ij 
(pnl) i p 
(Pnn)u 
{Pnl)ip 

(PLN)aJ 
(PLL)a/3 : 
(PLN)aJ 
(PLL)a/3 ' 



-■ (4,i a ^,i)/^ 

(a\ tl a/3 t i)/V 

■ (4,2 a JA)/ V 
(a\ t2 a/3 t i)/V 

(<&,iap,i)/V 



(Pnn)u 
(PnL)ip ■ 
(Pnn)u 
(PnL)ip ■ 

(PLN)aJ 
(PLl)aP z 
{PLn)»J 

(PLL)ap z 



■ (a\ i2 aj t2 )/V 

(4,2 a /3,2>/^ 
<<iOJ,2>/V 

(4,1«£,2>/^ 
( a t,2 a J,l)/ V 

( a L«/3,2>/^ 



(53) 



where we have suppressed time and momentum indices (all momenta are p and all times t). V is 
the overall spacial volume, which will always drop out of the computations in the end. 



29 Whether the ladder operators or p or both change with time depends on whether one choses the Heisenberg, 
Schrodinger or interaction picture. The expectation value, however, always has the same time dependence. 
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A. 2 Effective Kinetic Equations 

The time evolution of p is governed by an effective Hamiltonian. In absence of Hubble expansion, 
which we will add later, it follows the kinetic equation 

if t =[K,p]- l 1 {r > , P } + i{T<,l-p}. (54) 

H can be viewed as the dispersive part of the effective Hamiltonian. The absorbtive part given 
by the matrices arises because the system is coupled to the background plasma formed by all 
other degrees of freedom of the SM. Note that (|54p is valid for each momentum mode separately. 
The different modes are coupled by H and r<, which in principle depend on p and the lepton 
chemical potentials. The smallness of the sterile neutrino couplings F allows to simplify (|54p due 
to a separation of time scales: The time scale associated with the Nj dynamics and the time scale 
on which chemical equilibration of the lepton asymmetries occurs are much longer than the typical 
relaxation time to kinetic equilibrium in the SM plasma. This allows to employ a relaxation time 
appr oxim ation and relate T > and T < by a detailed balance (or Kubo-Martin-Schwinger) relation 

i^ = [B,p]- l -{T,p-p^}, (55) 

with r = r > — r < . p eq is p evaluated with an equilibrium density matrix, p = Z/trZ, Z = 
exp(—H/T), where H is the Hamiltonian corresponding to ([1]). The matrices H and T are 
Hermitian. 

The effective masses of active and sterile neutrinos are very different and fast oscillations 
between them play no role. We thus put to zero p NL , p LN , p^ L , p LN , p NL , p L ^, pf, L and p L ^. 
The time evolution of the asymmetries is related to the relaxation time scales of the TV/. Since 
interactions of the active neutrinos amongst themselves and with other SM fields are much faster, 
coherent effects in the active sector are negligible on this time scale. This allows to furthermore 
neglect p L i and pi L . pll and pn are taken diagonal with equilibrium occupation numbers and 
are thus characterized by the temperature T and three slowly varying chemical potentials. Thus, 
we can entirely describe the active sector by four numbers. Instead of the chemical potentials, 
we will in the following use n a = (pLi)aa — (pt r.)nn, i-e. the number of particles minus number 
of antiparticles, to characterize the asymmetrieqfj. The relation between both can be found in 



the appendix of 19(] . In the sterile sector we have to keep track of coherences. The system can 



then be described by the following set of kinetic equations, 

= 1 H ' PNNI ~ ^{1N,PNN ~ PnnS + 2 n <*lN i ( 56 J 

i dPM_ = ^psri-llfapw-pfay-ln^ , (57) 

= _ h2 n a + iTr[^(p NN -p^ N ))-iTr[jr(Pm-p e i^} ■ ( 58 ) 



We work in the F = FUn base in flavor space. This is the mass base of sterile neutrinos in vacuum, but the 
flavor base for active neutrinos. Hence, the diagonal elements of pnn and Pjvjv have a straightforward interpretation 
as number densities for physical particles in vacuum, while the ladder operators a' a create linear combinations of 
physical particles, and the matrices pll and pn are strictly speaking not diagonal in thermal equilibrium. This 
adds a subtlety to the interpretation of n a as "particles minus antiparticles", which, however, is of no practical 
relevance due to the smallness of the neutrino masses. 
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Here jn, 7^- are the appropriate block-diagonal submatrices of T, for the corresponding submatrix 
of H we used the same symbol as for the full matrix to simplify the notations. These equations 
do not take into account the expansion of the universe. As usual, it can be included by using 
abundances (or "yields") instead of number densities. It is also convenient to introduce the 
variable X = M/T rather than time t. 

All quantities appearing in the above equations depend on momentum. The different momen- 
tum modes are coupled by the scattering and decay processes. We have suppressed this momen- 
tum dependence. We define the abundances Pn = f d 3 p/(2Tr) 3 pnn / 's, pn = f d 3 p/ (2n) 3 p NN / s, 
pF9 = f d 3 p/(27r) 3 p eq N /s « / d 3 p/(27r fp e ^/s and y, a = / d 3 p/{2Ti) 3 n a /B Assumption 5) is 
justified if the common kinetic equilibrium assumption 



(pnn)ij (pn) 



1 j 



(P^ N )U (P eq ) 



1 j 



(59) 



holds. We can use ([59]) to rewrite the anticommutator in (f56|) as {Tn,Pn — p eq } with 



T N = T f<l 3 P7^=t f<i 3 P7^, (60) 
J P eq J n F 

n F = J <1VfK). (61) 

We again emphasize that pnn, In etc. appearing in ([56]) - ([58]) depend on momentum while Tj^, 
Pn, p eq etc. do not. 

For T <C M, almost all particles have the momentum p ~ T and is essentially obtained by 
evaluating 77V at p = p. Practically we compute the rates as described in section [3. 2. 2i Similarly, 
we can use H = x H| p= p for the Hermitian part of the effective Hamiltonian. For |p| ~ T > M, 
n F can be approximated by np |("(3)T 3 m 1.8T 3 , but has to be computed numerically. 

Using the above considerations, we can write down the following effective kinetic equations: 

= [H, PN ]- l -{T N ,p N -p^}+ l -p a t N , (62) 

4r = [ H *iPN\-\pN,PN-p eq }-\^ , (63) 

*W = - iT L^ + iTr[rUPN-p eq ))-iTr[rr(PN-p e9 )] ■ (64) 

They are equivalent to the ones used in [jj. Their interpretation is straightforward. In the 
mass base, the diagonal elements of pn and p N are the abundances of sterile neutrinos and 
antineutrinos, respectively. The off-diagonal elements are flavor coherences, p^v thus gives the 
abundances for "particles" and p N those for "antiparticles" , defined as the helicity states of the 
Majorana fields Nj. This interpretation holds in vacuum, while at finite temperature the effective 
mass matrix rotates due to the interplay between the (temperature dependent) Higgs expectation 
value, the Majorana mass Mm and thermal masses in the plasma. 

The first two terms in (j62[) and (j63|) are due to sterile neutrino oscillations and dissipative 
effects, respectively, either by scatterings or by decays and inverse decays of sterile neutrinos. 
More precisely, the Hermitian 2x2 matrix H in (|62p and (|63p is the dispersive part of the 



31 Note that the fi a defined this way are dimensionless and basically abundances, not chemical potentials, cf. 
appendix [C] 
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effective Hamiltonian for p^ and pn- The matrix T^v is the dissipative part of the effective 
Hamiltonian for p^ and pN that arises because the sterile neutrinos are coupled to the SM. p eq 
is the common equilibrium value of om and pat in absence of an asymmetry. All these terms 



appeared already in earlier studies 43]. The equations of motion f)64f) for the asymmetries in 
the active sector follow from consistency consideration and the symmetries of the i/MSM. The 
terms containing in (164p are their c ounte rparts in the active sector. The last term is due to 



backreaction and has been discussed in 1391 ] 



A. 3 The Effective Hamiltonian 

We follow the approach used in [8[] and split the Hamiltonian in the Heisenberg picture into a 
free part Hq and interaction Hint- We perform the computation in Minkowski spacetime and for 
the moment omit the factor dt/dX included in the definition (I60p . The same rates, multiplied by 
this factor, can be used in the early universe when abundances are considered instead of number 
densities. Starting point of the computation is the von Neumann equation in the interaction 
picture, 

where p\ = ex.p(iHot)pexp(—iHot) is the density matrix in the interaction picture and Hi = 
ex.p(iHQt)Hi nt exp(—iH()t), where p is the (time independent) density matrix in the Heisenberg 
picture. Equation ([65]) can be solved perturbatively, 

ft ft ft' 

h(t) = h-i / dt'[H I (t'),p ] + (-i) 2 / dt' / dt"[H 1 (t'),[H 1 (t"),p ]} + ..., (66) 
Jo Jo Jo 

where p = p{0) = p\(0). 

We use ([66]) to compute the expectation values (a\ r (p, t)aj >s (p, t))/V by insertion into ([28]) . 
For po we chose a product density matrix po = pj<i <8> p e gMi wnere Psm ls an equilibrium density 
matrix for the SM fields and 

p N = ^P/, s 4jO)(0|a /iS . (67) 
/ 

This is not the most general density matrix that can be build from one-particle Nj states, but it is 
sufficient to derive the effective Hamiltonian. The formula (]66p formally gives expressions for the 
bilinears at all times. These are strictly valid only at times much shorter than the sterile neutrino 
relaxation time because the perturbative expansion at some point breaks down due to secular 
terms. In the relaxation time approximation (j55h we can use a trick to deal with this problem. 
We differentiate (oj r (p, t)aj )S (p, t))/V with respect to time to obtain a "rate". We then send t 
to infinity to eliminate its explicit appearance from the rate. This last step is allowed because all 
correlation functions of SM fields are damped on time scales much shorter than the sterile neutrino 
relaxation time by thermal damping rates due to the gauge interactions. Thus, the late time part 
of the integrand in ([66]) does not contribute significantly to d(a\ r (p, t)aj jS (p, t))/(Vdt). This way 
we obtain the rate of change of the matrices pnn and Pfifj at initial time. In the relaxation time 
approximation, these rates can also be used at later times because backreaction is accounted for 
in the pn — P eq term. 
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Repeating literally all steps in section 2.2 of reference [8j for the two flavor case and the initial 
density matrix (|67|) . we obtain 



d (4,r Q J,s) 

Tt V ~ P> A 

A=l,2 



M 

i{$Ai ~ S AJ )—Re(AM M )ij 



1 l f J p iUJp) [(5ai + ^j)tr (P*(P- hjPLUjjip) + J^^j/ftllJ^)) 

+^A/-<W)P / ^^—tr{P R (P: s ) I jP L UJ J ( qo ,p) + P R (P: r )j I P L Uj I (q ,p)) 
J z-ir w p — f/ 



(68) 



with 



and 



AM M = M M - 1M , M M = UjjM M U N (69) 



(PDjj = u/, r (p)fij,,(p), (PD/j = uj,r(p)uj, a (p). (70) 
In the limit — > the projectors are independent of the sterile flavor indices and reduce to 

(P u )fj = (/+M)Q-(-l)V), (P v )fj = [J- M) Q + (-!)*{,) , (71) 

where f) is the helicity matrix (|5T1) . We have used u c = Cu T = v and introduced the self energies 

IT|» = F a iFpj I -0^ (v 2 5(p -k) + A^p + k)) S%(k), (72) 

with F = FUm- The thermal Wightman functions appearing therein are defined as 

A > (x 1 -x 2 ) = (<t>(xi)<t>(x 2 )) (73) 
A<(a;i-x 2 ) = (0(s 2 )0(si)) (74) 
S^pij{xi-X2) = (u ai (xi)u^(x 2 )) (75) 

S aPij( x l ~ X 2) = -(%(^2)^ai(^l))- (76) 

Here j,,- are spinor indices, which we suppress in the following. Transitions with r / s do not 
contribute at leading order in a j due to the projectors. This justifies our description of the sterile 
sector by two 2x2 matrices pjy and pfj rather than a 4 x 4 matrix including elements oc p^vjv 
etc. Transitions with a ^ ft are suppressed by the small active neutrino masses mi/T <C 1. The 
above expressions are written in the P-base (vacuum mass base). They can be translated into 
the F-base used in (UJ by the replacements F — >■ F and AMm — > AMa^. Note that M is defined 
at T = 0. 

With ([67|) . the initial value for p^ can be written as p^ oc diag(Pi, P 2 ). The RHS of ([68]) has 
a real and an imaginary part. They allow to extract the dispersive and dissipative parts H and 
7at of the effective Hamiltonian. 
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Figure 17: The Higgs expectation value as a function of temperature for mu = 126 GeV. 




A. 3.1 Dispersive Part H 

Comparison of (|68|) and (156p in absence of active lepton asymmetries (since we chose p £ g M without 
chemical potentials) allows to define the dispersive part of the effective Hamiltonian appearing 



m 



Hu 



T(M 2 +p 2 )^ /J + T 

np 



M 



Re(AM M )u 



(77) 



1 - 2f F (oj p ) 

2Wn 



V 



7n 



2ir uj t 



tr (p n n 7J (g , p) + P„II J7 (go, p)) 



where we have introduced the short notation 



P u = P R (Pl l )uP L 



(3 + *>)Pl 



V v = P r (P?)jiP l 



(78) 



The additional factor fp(ujp) /tif and the momentum integral come from the momentum averag- 
ing, cf (|60p . One can distinguish between three contributions. The term involving AMm comes 
from the splitting of the Majorana masses and remains present in vacuum. The term involving 



Hjj is due to the Yukawa interactions. It contains two contributions, see (|72p . which are related 
to the Feynman diagrams shown in Figure HI The part oc v(T) 2 is due to the interaction with 
the Higgs condensate and produces the Dirac mass at T < Tpw. The part involving comes 
from scatterings with Higgs particles. The Higgs expectation value as a function of temperature 
can b e ca l culat ed for a given Higgs mass. We used mrj = 126 GeV, as suggested by recent LHC 
data fl35l.fl36l]. to obtain the dependence shown in figure [TTJ However, we checked that varying 
mu within the allowed window 115 — 130 GeV does not have a big effect on the results. 

Evaluation of (|77p requires knowledge of the dressed active neutrino and Higgs propagators 
S<g(p) and A^, respectively. These are in principle complicated functions of p and T. However, 
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H T j = x(M 2 + T 2 )5 U + x J d 3 p- 



we are mostly interested in very high or low temperatures, T > Tew ^ M during baryogenesis 
and T < M in the context of DM production. This allows to simplify the expressions. It is 
convenient to dissect the self energy IT - into the Lorentz components 

P L U u (p) = P L (Au(pW+ BufrM, (79) 

where u = (1, 0, 0, 0) is the four-velocity of the primordial plasma, and write 

MwpH M ~ l-2/H^ p ) , . 

' — Re(AM M )u-\ — (80) 

rip ui p Zojp 

x \<p I ^£ _L_ (u p (B u + Bjj) + p(B u - Bn) + M 2 (A U + Ajij) 

Here the momentum dependence of Bjj(q Q ,p), Ajj(qo,p) has been suppressed and p has to be 
read as |p|. At temperatures T 3> M, the integral is dominated by hard momenta ~ T and the 
term involving Bjj dominates Hjj. For T < v(T), the interaction with the Higgs condensate 
dominates the Nj self energy and Bjj can be estimated as 

Bjj(p) ~ v(T) 2 F aI F*j ^-b af3 (p) (S(oj - |p| + b) - S(uj + |p| + b)) for T < v(T). (81) 



Here b is the so-called "potential contribution" to the active neutrino propagator [6|] , obtained by 
decomposing the retarded active neutrino self energy as 

ReSjaCp) = a aP {p)^+ b a p(p)yt. (82) 

Since active neutrinos mainly scatter via weak gauge interactions, the coefficients are in good 
approximation flavor independent in the primordial plasma, we can defi ne b5 a p = b a p(p), where 
b a s(p) is to be evaluated on-shell. For hard momenta, b gives 



8p \^ ~ cos 2 e w 
lggl . 9 „ \ 7tt 2 



2 + 7^7 , T->M W 



2 + cos 2 %) , T^M W . 



(83) 



Here is the Weinberg angle and aw the weak gauge coupling. This leads to a contribution to 
Hjj of the form 

v(T) 2 F aI F* 



lJ 2p M 2 + 26w p ' 



For T > Tevj/, scatterings with Higgs bosons dominate and the hard thermal loop result H ~ 
Ft FT/8 from |3j can be used. Combining the two contributions, we approximate H by 

H ~ —AM + (FF)* (| + ^) (84) 

during the calculations in section [5j In practice it is more convenient to work in the F-base, 
where the Hamiltonian reads -a 3 ^AM + F f F(| + ^Ip-). 

At temperatures T <C Tew , there are no Higgs particles in the plasma and the Higgs expecta- 
tion value is constant, thus Bjj = v 2 F*jFpjb a p, Ajj = v 2 F^ I Fpja a p. In [f| it has been estimated 
that thermal corrections to the active neutrino propagator are small below a temperature 



M \ 3 

Tpo * = 13 lGev) GeV ' (85) 



(. 
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For the masses under consideration in this work, we can approximately use free active neutrino 
propagators in section [6] because T + < T pot . Furthermore, due to the considerations in section 
2.61 we are mainly interested in the case Un ~ 1 for DM production, thus F ~ F. Then 
Wpo,|p|) =^0and a a p(p ,\p\) ~ (U v )ai(U v )*^(S(po - Ui) - 6(p with = (p 2 +mf) 1 / 2 



1421 ] . where mj are the active neutrino masses. This recovers the vacuum result for the mass 



matrix at |p| = 0, cf ([9]). H can be approximated by 

H = x(p 2 + M 2 N ) X I 2 . (86) 

In the basis of vacuum mass eigenstates it has the form H = diag((p 2 + M|) 1//2 , (p 2 + M 2 ) 1 / 2 ). 
Since 5M <CT,Mwe can expand in 5M and obtain for p = T 

H — x 5M(X~ 2 + l)~ 1/2 a 3 , (87) 

with X = M/T and the third Pauli matrix 173. The part of H that is proportional to the identity 
matrix has been dropped as it always cancels out of the commutators in the kinetic equations. 

A. 3. 2 Dissipative Part 

Again comparing (j68[) and (|56p . we define 

(Fs)lJ = T / rf 3 p /i?K) 1 = 2 / FK) tr (P u U u (p) + PJIJjfr)) ■ 
J TIF 

The rate for the "antiparticles" can be found by using projectors analogue to ([78]) . but with 



helicity index 2- It is given by (r^)* as expected 32 !. For what follows, it is useful to pull the 
Yukawa matrices out of the self energies and define 

D§j(p)= / ^(v 2 6(p-k) + A^(p + k))S%(k). (89) 



Obviously 11^ = F^jFpjU^p. For the computation of Tat according to ([39]) we can now define 
the matrices 

R(T,MU = J <* 3 P /F ^ P) 1 " 2 2 ^ K) tr (P^fr)) (90) 

Rm(T, M) a p = J 1 ~ ^^ tr (P„fi^(p)) ■ (91) 

They in general have to be computed numerically. We discuss their properties in section 14.21 
As usual in thermal field theory, the sterile neutrino self energies II < and II > can be associated 



N 



with the gain and loss rate. Their difference IT - = II > — II < gives the total relaxation rate T 
for the sterile neutrinos. It acts as thermal production rate when their occupation numbers are 
below their equilibrium values and as dissipation rate in the opposite case. In configuration space, 
the self energy II - (:z;) is related to the retarded self energy by U R (x) = 9(xq)H~ (x). This implies 
n~(p) = 2ilmTl R (p) in momentum space. As usual in field theory, the imaginary part of YL R can 
be related to the total scattering cross section by the optical theorem (or its finite temperature 



32 This can be seen by noticing that the traces are real and PrP^ 2 Pl = PrPI Pl, PrPv 2 Pl = PrPu 1 Pl under 
the trace. 
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generalization) 1281 ] . while the real part is responsible for the mass shift (or modified dispersion 



relation in the plasma). Both are related by the Kramers-Kronig relations. The appearance of 
II in ([88]) is in accord with the optical theorem, and the contributions to the dispersive and 
dissipative parts of the effective Hamiltonian are indeed related by a Kramers-Kronig relation, cf 
([77]) and ([88]) . This provides a good cross-check for our result. 

A. 3. 3 The remaining Rates 



The remaining rates (]40p and (j4ip appearing in (|29j) - f)31 j) in principle have to be calculated 
independently. The precise computation is considerably more involved than in the case of T^. 
Tjv is related to the discontinuities of the Nj self energies, which to leading order in the tiny 
Yukawa couplings F a j only contain propagators of SM- fields as internal lines. Due to the fast 
gauge interactions these are in equilibrium in the relaxation time approximation and the RHS 
of (|88p can be computed by means of thermal (equilibrium) field theory. This is not possible in 
the computation of the damping rates for the SM-lepton asymmetries, which are related to self 
energies where the out of equilibrium fields Nj appear as internal lines. For simplicity we follow 
the approach taken in 0], see section 6 therein, and use the symmetries of the i^MSM in certain 
limits to fix the structure of the rates. 

We first consider the limit (Mm)ij —> 0, that is the absence of a Majorana mass term. In 
this case the "total lepton number" is conserved, 



\ I a ) 



o = %\y A + y (92) 

M M =0 

4 = VRffviU (93) 
J% = VL, a ^VL, a + e.L, a l^e L , a . (94) 

To leading order in the small mixing 8 a j this implies 




, fl>a ) -0. (95) 

M M =0 

This situation is in good approximation realized for T S> M, when baryogenesis takes place - the 
total lepton number is not violated during this process and a non-zero baryon number is only 
realized because sphalerons couple exclusively to left handed fields. Equation ([95]) implies 



r^v-Eafg trf^~T2 for(M M )/j = 0. (96) 



Other interesting limits considered in [6J include F a j — > with fixed a for all I (leading to 
conservation of J a) and F a j — > with fixed I for all a (leading to individual conservation of the 
combination Jj + Ja and the remaining current Jj,j). These limits allow to fix the basic 
structure of equations (]40p . (|4ip . For a general choice of parameters some corrections of 0[1] to 
these relations may be necessary, the determination of which we postpone until the precision of 
experimental data on the z^MSM requires it. 

A. 4 Uncertainties 

Our study is the most complete quantitative study of bounds on the z^MSM parameter space 
from cosmology to date. However, the various assumptions made in the derivation of the kinetic 
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equations lead to uncertainties that may be of order one. These can be grouped into three 
categories: 



We only consider momentum averaged quantities. Since the sterile neutrinos can be far 
from thermal equilibrium, one in principle has to study the time evolution of each mode 
separately. Our treatment is a reasonable approximation as l ong as the kinetic equilibrium 
assumption ([59]) holds. A study of this aspect published in [1391 ] suggests that deviations 
from kinetic equilibrium are indeed only of order one . 



The rates (j40|) and (|4ip have been calculated in a rather crude way in section IA.3.31 leading 
to another source of uncertainties of order one. In addition, a precise calculation of the BAU 
requires knowledge of the sphaleron rate throughout the ele ctrow eak transition. Including 
this is expected to yield a slightly bigger value for the BAU 1431 ] . 



• Though they are matrix valued and allow to study flavor oscillations, the equations ([62p ~(j64 
are of the Boltzmann type. They assume that the system can be described as a collection of 
(possibly entangled) individual particles that move freely between isolated scatterings and 
carry essentially no knowledge about previous interactions ("molecular chaos"). 

The first two issues can be fixed by more precise computations. However, with the current exper- 
imental data, order one uncertainties are small compared to the experimental and observational 
bounds on the model parameters. The corrections will only slightly change the boundaries of 
the allowed regions in parameter space found in this work. We therefore postpone more precise 
calculation to the time when such precision is required from the experimental side. 

In contrast to that, the third point is more conceptual. In a dense plasma, multiple scatterings, 
off-shell and memory effects may affect the dynamics. The effect of these cannot be estimated 
within the framework of Boltzmann type equations, it requires a derivation from first principles 
that either confirms f[62|) - ()64[) and allows to estimate the size of the corrections or replaces them 
by a modified set of equations. In the past years, much pr ogress has been made in the d erivation 
of effective kinetic equations from first principles 75-78, 125 . 1261 . 130 . 1331 . 144 - 168 ]. Recent 
studies suggest that kinetic e quation of the Boltzmann type are in principle applicable to study 
leptogenesis 77, [71) 154 . 168 ]. but the resonant amplification may be weaker than found in the 
standard Boltzmann approach 77j]. It remains to be seen which effects possible corrections have 
in the i^MSM, where baryogenesis and dark matter production both crucially rely on the resonant 
amplification. A first principles study is difficult in the i/MSM due to the various different time 
dependent scales related to production, oscillations, freezeout and decay of the sterile neutrinos 
and the vast range of relevant temperatures. However, at this stage it seems likely that a first 
principles treatment is, if at all, only of phenomenological interest in the region around Tjv ~ 5M, 
which makes up only a small fraction of the relevant parameter space, cf. figure [H 



B Connection to Pseudo-Dirac Base 

In our notation, the elements of pat and are bilinears in ladder operators that create quanta 
of the fields Nj, i.e. mass eigenstates in vacuum. This has the advantage that the diagonal 
elements can be interpreted as abundances of physical particles. The rates R and Rm have been 
introduced in 0], to which we regularly refer in this article. The basis in the field space right 
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handed neutrinos there differs from the one we use in (p} and corresponds to Uvr with 

In this basis Mm is not diagonal and the Yukawa couplings should be rewritten as hU = F. The 
computation of the rates is then performed by defining a Dirac-spinor 

* = U21VRJ + (U 3I u RJ ) c . (98) 

This is possible when the small mass splitting between the sterile neutrinos is neglected (or viewed 
as a perturbation and placed in the interacting part of C). The fields vrj can be recovered from 
this as vrj = U^Pr^ + U 3I Pr^ c . In terms of the i/MSM Lagrangian reads 

C = £ S M+£o + £int (99) 

-iAAf(W + W c ). (100) 
The analogue of our matrix p?q (which is also called pj^ in 0]) is defined as 

_ f d 3 p ( (c\ Cl ) <cK> \ 

PN ~ J w ^ (4ci) (d\ dl ) ) ' (101) 

where c s (cj) and d s (d\) are the annihilation (creation) operators for ^ particles (antiparticles) 
with momentum p and helicity s. The corresponding rate in the kinetic equations is given 
bj@ 



1 N 



t J d 3 p^^v 2 ^^(h aI h}jR(T,MU 

where h = FU\ a± is the first Pauli matrix and we have neglected flavor off-diagonal elements. 
In the high temperature regime that was considered in 0] this simplifies to 

T^t" 1 = (rfh)*R( s XT,M) + aitfhaiR { M > (T,M). (103) 

a±h^ha\ is (h^h)* with the diagonal elements swapped. The equation (j90|) defines the quantities 
R(T, M) and Rm(T, Aff^l. Ignoring the small mixing between active and sterile neutrinos, T^y is 
related to by 

r N ~ (uu N ) T r%(uu N y (104) 

for T 3> rrii. 

Finally, the Yukawa matrix in and 0] is expressed in terms of the parameters e du ri di , <f> di , 
which differ from those we use here. In the limit e d i >> 1, these can be related to our parameters 
by 

y/eJi = e~ lmw , r] di = 2Recu, a di = 4> di = 5. (105) 

We here prefer to use the parameterization fixed in section [231 because the expressions in therms 
of (|105p given in 0] are only approximate. 



33 Note that there were errors in the corresponding expressions (5.19), (5.20) in [g] and that there only the case 
T> If was considered. 



34 Our definitions of R and Rm differ from those in [f| by a constant factor F§ with Fq = 2 x 10" 



!) 
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C How to characterize the lepton asymmetries 



In the z^MSM neither the individual lepton numbers, related to the currents (I94p . nor their sum 
are conserved. However, since the rates of all processes that violate them are suppressed by the 
small Yukawa couplings F, they evolve on a much slower time scale than other processes in the 
primordial plasma and are well-defined. For practical purposes the magnitude of flavoured lepton 
asymmetries in the primordial plasma can be characterized in different ways. In this article we 
describe them by the ratio between the number densities (particles minus antiparticles) and the 
entropy density s = 2TT 2 T 3 g*/ '45, 

li a = — cf. d. 

s 

This quantity is convenient because it is not affected by the expansion of the universe as long as 
the expansion is adiabatic. In the following we relate [i a to other quantities that are commonly 
used in the literature, using the relations given in the appendix of [191 ] . 

In quantum field theory calculations it is common to parameterize the asymmetries by chem- 
ical potentials |J. a , which can be extracted from the distribution functions that appear in the free 
propagators at finite temperature. In the massless limit T 3> m; these are related to n a by 



u T 2 u 3 

n a = — h— -o, (106) 

O 07T Z 



leading to 

^A^T- (107) 

Alternatively one can normalize the lepton numbers n a ("particles minus antiparticles") by 
the total density of "particles plus untiparticles" in the plasma, 

A a = ^ (108) 
where n a q = 2 / d 3 q/(2^) 3 /(e |q|/T + 1) = 3C(3)T 3 /2^ 2 and 

,« = ^|fW (109) 

Finally, one can normalize with respect to the photon density, 

L a = ^ (110) 

n 7 

where n 7 = 2 / d 3 q/(2vr) 3 /(el q l/ T - 1) = 2C(3)T 3 /vr 2 , which yields 



,„ = ^«. an) 



D Low temperature Decay Rates for sterile Neutrinos 

Most of the rates relevant for this work have been computed in where they are listed in the 
appendix. Here we only list those rates that are needed in addition to those or require refinement. 
This was necessary for the decay rates into leptons, where masses of the final state particles had 
been neglected in the original computation. 
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D.l Semileptonic decay 

Decay into up-type quarks through neutral current : 



Nx^fV a UU 



G 2 F \e aI \ 2 M 5 (, u) 



1927T 3 



16 



fM(x q )S(x q ,x q )+x* ( 3 - — C x x\ + (3 - 8d) 



x log 



1 - 4x 2 q + 2x\ + S(x q ,x q )(l - 2x 2 q ) 



2x\ 



(112) 



where x q = m q /M. 

Decay into down-type quarks through neutral current : 



_ G 2 F \Q aI \ 2 M 5 

N^Vadd 1927T 3 



fW(x q )S(x q , x q )+x q ( 3 - \c 2 x* - (1 - ic 2 )x* 



x log 



l-Ax 2 q + 2x i q + S{x q ,x q )(l-2x 2 q ) 



2x q 



(113) 



Decay into quarks through charged current : 
_G 2 F \V nm \ 2 \e aI \ 2 M 5 ' 



iV/— >e a u n dr, 



1927T 3 



g(x,y)S(x,y) 



12x 4 log 



12y 4 log 



+ 12x 4 y 4 log 



1 - S(x, y)(l + x 2 - y 2 ) - 2y 2 + (x 2 - y 2 ^ 2 



2x 2 



1 - S(x, y)(l -x 2 + y 2 ) - 2x 2 + (x 2 - y 2 ) 2 ^ 



2y 2 



1 - 2x 2 - 2y 2 + x 4 + y 4 - S(x, y)(l - x 2 - y 2 ) 1 
2x 2 y 2 



(114) 



where min(m Q , m Un , md m ) is neglected, and x and y are the two heavier masses divided by M. 
D.2 Leptonic decay 



G 2 F \Q aI \ 2 M 5 



192vr 3 



12x 4 log 



-12xt log 



S(x a ,xp)g(x a ,xp) 



1 - S(x a ,xp)(l + x 2 a - xf) - 2x| + {x 2 a - x 2 p ) 2 
2^ 

1 - S(x a , Xp){\ - X j+ X 2 ) - 2x 2 a + (4 - tf}? 
2T 2 



+ 12x 4 x 4 ,log 



(115) 
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with 



S(x,y) 
Ci 



V(l-(x + y) 2 )(l-(x-y)2) 
sin 2 Ow (3 -4 sin 2 W ) 



I ~ \C, - Q - ^d) x 2 - Q + 4d) x 4 + (-3 + 8d) x 



sin 2 (3 — 2 sin 2 




11 in 2 i 

I - g C2 + " 7 )X2 " ( 2 + 2Ca)x4 " (3 " 4Ca)x6 

1 - 7x 2 - 7y 2 - 7x 4 - 7y A + 12x 2 y 2 - 7x 2 y 4 - 7x 4 y 2 + x 6 + y 6 . 
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